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ABSTRACT 

Specific  heat  measurements  on  three  samples  of  polyethylene, 
differing  only  in  density,  were  made  in  the  temperature  range  of  2.H 
to  30°K.   A  definite  density  dependence  was  noted  for  the  specific 
heat  in  this  temperature  interval  which  allowed  extrapolation  of  the 
data  to  completely  crystalline  and  completely  amorphous  cases.   At 
the  lowest  temperatures,  the  amorphous  results  were  observed  to  dis- 
play an  "excess"  heat  capacity  which  could  be  accounted  for  by  the 
occurance  of  a  single  delta  function  peak  in  the  low  frequency  part 
of  the  vibrational  spectrum.   This  excess  did  not  appear  in  the  com- 
pletely crystalline  extrapolated  data,  and  the  specific  heat  was  found 
to  be  proportional  to  the  cube  of  the  temperature  up  to  9°K. 

An  attempt  was  made  to  compare  the  results  with  two  theoretical 
models  (Tarasov  and  Stockmayer-Hecht)  with  only  fair  agreement  in  one 
case,  and  none  in  the  other.   The  agreement  with  previous  experimental 
results  is  excellent  if  the  density  dependence  and  the  excess  heat 
capacity  are  considered. 
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1 .   INTRODUCTION 

General  Remarks 
The  objective  of  the  investigation  which  is  the  subject  of  this 
thesis  was  to  determine  the  specific  heat  of  polyethylene  by  direct 
calorimetric  methods  in  the  temperature  range  between  2.4  and  30°K. 
This  study  was  prompted  by  some  earlier  results  of  low  precision  which 
indicated  a  rather  sizeable  density  dependence  of  the  specific  heat  in 
the  low  temperature  range  [34]  .   This  observation  was  in  contradiction 
to  statements  in  the  literature  which  denied  such  a  dependence  [52]  . 
In  addition,  theoretical  models  had  predicted  an  unusual  behavior  of 
the  specific  heat  of  fibrous  solids,  such  as  polyethylene,  at  low 
temperatures  which  had  not  been  verified  by  previous  measurements 
[l9,  46,  47,  52]  .   It  was  hoped  that  these  effects,  if  they  were  to 
occur,  would  be  found  in  the  range  covered  by  the  present  experiment 
thus  affording  a  comparison  of  the  experimental  results  with  the  pre- 
dictions.  Further,  it  was  desired  to  investigate  the  density  depend- 
ence over  a  larger  temperature  range,  and  with  more  precision  than 
that  in  the  work  noted  above. 

The  unusual  heat  capacity  behavior  which  had  been  predicted  to 
take  place  is  a  consequence  of  the  polyethylene  molecules  being  long 
chains,  and  the  forces  along  the  chains  being  much  larger  than  those 
between  chains.   This  one  dimensional  characteristic  gives  rise  to 
theoretical  descriptions  for  the  specific  heat  which  differs  from 
most  other  non-fibrous  solids.   The  theoretical  models  based  on  this 
chain  structure  will  be  discussed  later,  as  will  the  morphology  of 
polyethylene.   However,  for  now  it  should  be  noted  that  polyethylene 
is  an  example  of  a  polymer  which  crystallizes  partially,  and  thus  is 
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neither  completely  crystalline  nor  completely  amorphous.   So,  one 
parameter  in  any  discussion  of  polyethylene  must  be  of  its  crystal- 
linity,  which  is  directly  related  to  its  density.   Making  the  assump- 
tion that  polyethylene  consists  of  two  distinct  forms,  the  crystalline 
and  the  amorphous,  it  is  much  simpler,  for  this  discussion,  to  des- 
cribe the  specific  heat  of  polyethylene  in  terms  of  the  specific  heat 
of  these  two  idealized  phases,  i.e.,  the  specific  heat  of  any  real 
sample  can  be  expressed  as  an  admixture  of  that  of  the  two  hypotheti- 
cal ones.   In  order  to  accomplish  this,  the  heat  capacity  of  various 
samples  differing  in  crystallinity  (density)  was  determined  so  that 
extrapolation  to  the  completely  amorphous  and  crystalline  cases  could 
be  made. 

Heat  capacity  data  for  polyethylene  is  available  from  various 
sources  in  the  temperature  ranges  from  1  to  5°K  and  from  20°K  to  its 
melting  point  (about  428°K) .   The  results  of  this  research  will 
bridge  the  existing  temperature  gap  and  should  be  of  interest  for  two 
reasons.   First,  since  polyethylene  is  an  archtype  of  the  chain 
polymers,  the  behavior  of  its  heat  capacity  over  the  entire  tempera- 
ture might  be  representative  of  polymers  in  general,  and  second,  to 
allow  a  possible  check  on  the  aforementioned  theoretical  predictions. 

Previous  Measurements 

Below  100°K  there  is  not  much  experimental  data  on  the  specific 
heat  of  polyethylene.   Above  this  temperature  a  considerable  number 
of  measurements  have  been  made  to  determine  the  specific  heat  of  the 
solid  as  it  passes  through  the  glass  transition  and  the  melting  point. 
However  these  latter  results  are  not  germane  to  the  present  investiga- 
tion. 
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The  data  for  the  temperature  range  below  5°K  were  available  from 
two  sources.   Isaacs  and  Garland  reported  measurements  on  a  single 
high  density  sample  in  the  1.8-5. 3°K  interval  obtained  by  calorimetric 
methods  [21J.  The  heat  capacity  of  three  samples  of  various  densities 
was  measured  by  indirect  means  in  the  range  of  from  1  to  4.5°K  by 
Reese  and  Tucker  [34J  .   The  density  dependence  noted  before  is  quite 
apparent  when  all  four  sets  of  measurements  are  compared  [$4]  .      It 
should  be  noted,  however,  that  the  data  was  analyzed  assuming  a  cubic 
temperature  dependence  for  the  specific  heat.   This  assumption  is  not 
verified  by  the  present  investigation,  which  shows  a  distinct  non- 
cubic  relation  in  this  temperature  region;  however,  the  density  de- 
pendence is  valid. 

Dainton,  Hoare,  Evans,  and  Melia  give  results  for  two  samples  in 
the  range  20-300°K  [l3j  .   However  no  conclusions  could  be  made  about 
the  density  dependence,  since  the  samples  were  very  similar  in  density. 
The  only  other  known  source  of  specific  heat  data  for  polyethylene  in 
this  temperature  range  is  the  work  of  Sochava  from  17  to  60°K  [44 ]. 
The  sample  used  by  Sochava  was  not  well  characterized,  the  only  infor- 
mation given  being  that  it  was  a  "typical  commercial  product".   Thus 
there  is  no  way  to  determine  the  density  relationships  in  the  higher 
temperature  range  above  20°K.   Above  90°K  there  are  more  measure- 
ments but  it  is  felt  that  this  information  was  not  of  interest  insofar 
as  the  present  investigation  is  concerned. 

Brief  Description  of  Polyethylene 

Chemically,  polyethylene  is  one  of  the  simplest  forms  of  chain 
polymers,  consisting  only  of  CH2  units.   The  bonding  along  the  chain 
is  covalent.   The  chain  length  can  be  varied  during  manufacture,  but 
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in  general   it   is  quite   long,   on  the  order  of   1000  units.      In  addition, 
side  branches   can  and  do   form  on  the  chain,    and   is   also   a  parameter 
which  can  be  varied   according  to  the  manufacturing  process.      These 
branches,   which  are  short  chain   sections  normally  only  2-5   units 
long,    join  to  the  main  chain  by  the  substitution  of   a  carbon-carbon 
bond   for   a  carbon-hydrogen  bond    [31,38]  . 

Polyethylene   is  one  of  the  polymers  which  crystallizes  partially 
in  the  bulk  form.      A  necessary  condition  for  crystallization  is  geo- 
metrical  regularity  of  the  molecular  structure,   which   is   fulfilled  by 
the  simple  form  of  the  polyethylene  chain    [l@,3lj.      This  regularity 
allows  the  chains  to   form  ordered   arrays,   which  are  considered  the 
crystalline  state.      However  on  the  basis  of   single  crystals  grown 
from  solutions,   which   are  quite   small,    it   appears  that  the  crystals 
are  lamella,    approximately  100  Angstrom  units  thick,  with  the  molecu- 
lar chain  axis  being   perpendicular  to  the  face    [23, 35 J    .      whether  this 
crystalline  form  is  the  same   in  the  bulk   is   not   known  conclusively  [  31  ] 
However,    it   is  well  established   that  the  crystalline  part   of  poly- 
ethylene is  composed  of   clusters  of  spherulites   [ 31  J  . 

The  most   important  factor   inhibiting  crystal   formation   is  branch- 
ing.     This   is   evidenced  by  the   fact  that  highly  branched  polyethylene 
displays   relatively  low  crystalline  content  [_^+5 J   .      This   is   quite 
reasonable  since  the  regularity  of  the  structure   is  broken  at  the 
branch  points,   thereby  producing  discontinuities  or  distortions   in 
the  crystal. 

If  the  chain   is  not   in  the  well   ordered  array  of  the  crystalline 
state,    it   is   in  some  metastable  position,    intertwined  and  twisted 
around  other  chains.      The  complexity  of  the  non-crystalline  form  makes 
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description  very  difficult.  However,  it  sufficient  for  this  discus- 
sion to  characterize  the  non-crystalline  part  of  the  polyethylene  as 
the  amorphous  regions. 

Since  in  the  crystalline  part  of  the  polyethylene  the  molecules 
are  more  closely  packed  than  in  the  amorphous,    a  higher  density  is 
found   for  the  crystalline  phase   than  for  the   amorphous  [l5J  .      In 
fact   if  the  density,  Q    ,  of  a  particular  sample  is  known,    its  weight 
percent  crystallinity,  /   ,   can  be  calculated  from  the  relation, 


1.1 


where    />c    and    P/\  refer  to  the  densities  of  the  crystalline  and 
amorphous   parts  which  are  abtained  by   indirect  methods  f 15 J   .      The 
crystal linities  of   bulk  polyethylene  vary   in  general  between  35  and 
95%. 

Heat  Capacity  Models 
The  specific  heat  of  a  solid  can  be  calculated   if  the  tempera- 
ture dependence  of   its  internal   energy ,1^,    is  known.     This  follows  from 
the  definition 

1.2  Cv   "    P  T/v  ' 

If  the  internal  energy's  temperature  dependence  is  due  to  vibrations 
of  the  constituent   parts  of  the  solid  about   their  equilibrium  posi- 
tions,  the  total  energy  of  oscillation  over  the  whole  range  of  the 
vibrational   spectrum  is, 


1.3 


m=   f' *y(*)  Ef(*>r) Jv  , 
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where    £/»(^vr    "#  v  Av  /      is  the  average  energy  of  a  Planck 

oscillator  with  vibrational  frequency  "V,  and  °k  LV)  is  the  distribution 
function  for  the  eigenfrequencies  of  elastic  vibrations  of  the  solid, 
which  is  called  the  frequency  spectrum  [  14]  .  If  there  are  N  mass 
points  in  the  solid,  then  there  can  be  only  3N  normal  modes  of  vibra- 
tion, assuming  that  the  temperatures  are  sufficiently  low  that  possible 
internal  motion  of  the  mass  units  can  be  neglected  (the  optical  modes 
of  vibration).      This   sets  the  restriction 

Rigorous  treatment   shows  that  the  upper  limits   in  the   integrals   in 
equations   1.3  and   1.4  can  be  replaced  by  a  maximum  or  cut-off   fre- 
quency V^   !_50]'      Hence   if   g.  (y)   is  known,    "^    can  be  determined   from 
equation  1.4. 

This  cut-off  frequency  can  be   interpreted  as  the  frequency  at 
which  two   neighboring  mass  units  differ   in  phase  by   180°.      It   should 
also  be  noted  that   ")£    is  related  to  the  forces  acting  between  the 
neighbors,   that   is,   greater   forces   imply   larger  cut-off  frequencies, 
assuming  particles  of  equal  mass  [  14  ]  . 

The  central   problem  in  the  calculation  of  the   internal  energy 
of  a  system,   and  therefore   its   specific  heat,    is  to   find  the  frequency 
spectrum  <v^V)  •      This   is   in  general  extremely  difficult.      Therefore, 
several   approximation  methods  have  been  employed   for  solids.      One  of 
the  most   successful   approximations   is  the  Debye  model. 

The  Debye  approximations   are  that  the  frequency   spectrum  is  that 
calculated  assuming  the  solid  is   an  isotropic  continuum,    subject   only 
to   the  condition  stated  above  concerning  the  cut-off   frequency,    and 
that  there  is  a  linear  relation  between  the  frequency  and  the  wave 
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number,  i.e.,  that  there  is  no  dispersion.  These  approximations  can 
be  expected  to  hold  for  long  wavelengths  (low  temperature)  [  50  J  . 
Debye  theory  than  predicts  a  frequency  distribution  given  by 


1.5        <j(tt)Jv- 


pr*,  -)T  <y^ 


Vi  <\r  > 


Using  this  in  equations  1.2,  1.3  and  1.4  leads  to  the  familiar  Debye 
equation  for  the  specific  heat, 

i.6        <r„  -  3  R  P*  (¥/ 

where   P3  (±f)   »  / 2  ( ^  J#   -pCf      3   £"^-| 

and     vjj  -   -tt"     is  a  characteristic  temperature  of  the  material  called 

-r  i 

the  Debye  temperature.      For  very  low  temperatures,    -r^     ^  "^T"    , 
equation  1.6  reduces  to, 

1.7  C*  ~      5K        g-       y-&0/ 

and  for  very  high  temperatures, 

1.8     C^IR  T-»<*< 

This  last  result  is  the  same  as  is  predicted  classically  using  the 
equipartition  of  energy  [ 24  J  . 

As  a  logical  extention  of  this  model,  it  is  more  reasonable  to 
assume  that  there  exists  a  cutoff  wavelength,  rather  than  a  maximum 
frequency.   Since  it  is  known  that  for  all  solids  the  velocity  of 
sound  is  dependent  of  the  direction  of  propagation,  the  cut-off  fre- 
quency will  also  be  different  for  the  different  modes  and  directions 
L 43  J.   Recognizing  this,  equations  1.5  and  1.6  become, 
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3h/  (^f  ~vfl)aiv~  v-<^ 


w        y* 


iff  ^^nF  v;<^<^ 


1.10 


C,/*    fc 


M^+M*)] 


where  the  subscripts  refer  to  the  transverse  and  longitudinal  modes, 
and  the  y^s     are  directly  related  to  the  average  sound  velocities  in 
these  two  modes.   It  should  be  noted  that  this  modification  preserves 
the  T^  dependence  of  the  specific  heat  at  the  lowest  temperatures. 
Polymers  differ  from  most  other  solids  in  the  sense  that  they 
are  essentially  one  dimensional  in  nature,  because  the  forces  acting 
along  the  chains  are  much  greater  than  those  between  chains.   If  a 
Debye-type  analysis  is  made  for  a  one  dimensional  system,  in  a  manner 
completely  analogous  to  the  preceeding,  one  obtains  the  following  re- 
sults [  47]; 


where 


and     -&,   =  4j£  , 

The  model  employed  above  describes  non-interacting  chains,  so 
is  not  of  direct  interest  in  this  discussion.   However,  if  the  chains 
are  intertwined  and  linked  with  one  another,  as  in  a  chain  polymer 
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such  as  polyethylene,  there  will  be  both  the  three  dimensional  and  one 
dimensional  continua  existing  together.  This  is  the  basis  of  a  model 
due  to  Tarasov[47]  .   Since  the  forces  acting  along  the  chain  are  much 
greater  than  those  between  chains,  it  is  expected  that  the  maximum 
frequency  for  the  one  dimensional  part  of  the  model  will  be  greater 
than  that  for  the  three  dimensional  part.  Tarasov  therefore  predicted 
that  there  would  exist  two  frequency  spectra  such  that, 


1.13 


3 A/-  J  c^l(r)dv-    -M  <%,(vl& 


'^ 


where  the  subscripts  refer  to  the  one  and  three  dimensional  quantities. 
Using  the  Debye  assumptions,  frequency  distributions  can  be  found. 
However,  the  separation  into  transverse  and  longitudinal  parts  is 
again  applicable,  as  it  was  above.   Fbr  the  remainder  of  this  deriva- 
tion only  the  longitudinal  part  will  be  considered,  and  the  transverse 
contributions  put  in  at  the  end  by  analogy  with  the  Debye  case.   There- 
fore equation  1.13  becomes, 


1.14 


for  the  longitudinal  mode  only   (denoted  by  the  superscript  J?  ).     This 
equation  can  be  rewritten  as 

1.15  N    =    ^36     +  V,3 

where  N30   and  N^3  refer  to  the   first   and  second   integrals   above. 
The  distribution  of  the  modes   are,    for  N^3 


1.16  f 


1TJ<  )T  •     \r<YfJ 
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and,  for  N30, 


ij»  v  <  V* 


o  vAv. 


If  these  spectra  are  the  correct  ones,  the  heat  capacity  is  given  by 
the  resulting  Tarasov  equation, 

,18  cy-^(¥)-^[o^)-^m) 

for  the  longitudinal  mode,  where  "wu  {®Hl   are  related  to  V~,    j  Vj 
and  the  other  symbols  have  been  previously  defined.   Inclusion  of 
the  transverse  modes  then  yields  the  complete  Tarasov  relation, 

1.19  cv  =  c/  +  a(?J  . 

In  the  two  temperature  extremes,    equation   1.19  reduces   to 

1.21  Cv-    Jn  ' 

Between  these  extremes,    there  will   be  a  region  for  which, 


1.22 


Cv=  3i< 


No  continuum  model  should  be  expected  to  explain  the  specific 
heat  in  the  presence  of  dispersion,  which  will  take  place  in  any  real 
lattice.   In  order  to  calculate  the  density  of  states  in  the  vibra- 
tional spectrum  correctly,  the  effects  of  all  of  the  forces  acting 
between  the  mass  units  must  be  included.   In  practice  this  is  not 
feasible.   As  a  result,  the  usual  method  of  approach  is  to  make  real- 
istic estimates  of  the  strength  of  the  forces,  and  drop  all  but  the 
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dominate  ones  from  the  calculations.   Blackman,  and  many  others  since, 
have  achieved  a  certain  amount  of  success  along  these  lines;  the  de- 
gree of  success  depending  on  the  validity  of  the  assumptions  about 
the  relative  strength  of  the  forces,  and  the  number  of  forces  included 
in  the  calculations  [  6,7,16,28].   This  type  of  analysis,  since  it  em- 
ploys the  discrete  nature  of  the  solid,  should  be  more  complete  than 
the  continuum  treatment,  if  the  correct  parameters  are  used.   There  is 
however  one  feature  of  the  lattice  dynamics  calculations  which  is  the 
same  as  that  obtained  by  the  continuum  assumptions:  the  heat  capacity 
at  the  lowest  temperatures  should  be  proportional  to  the  cube  of  the 
absolute  temperature. 

For  crystalline  chain  polymers  an  attempt  has  been  made  to  obtain 
the  vibrational  spectrum  using  lattice  dynamics  by  Stockmayer  and 
Hecht  whose  work  was  greatly  extended  by  Genensky  and  Newell  fl9,46J. 
This  analysis  is  based  on  the  consideration  that  there  are  strong 
forces  along  the  chains  and  weaker  ones  between  them.   The  frequency 
spectrum  deduced  was  quite  complicated,  but  has  certain  characteris- 
tics which  are  predicted  to  be  characteristic  of  the  heat  capacity  of 
fibrous  solids  in  the  crystalline  form,  and  thus  should  be  applicable 
to  crystalline  polyethylene.   These  predictions  include,  (where  C^ 
now  refers  to  the  specific  heat  contribution  from  the  interchain 
interaction,  and  C3  to  that  along  the  chain)  [l9] 

ac.  +  c,. 

1.23  Cv  =      3 
and 

1.24  C,  «    T"*  -T«  ^  Xm 

1.25  C,  oc  T-'*  -^Xh«T«^Tw 


1.26  C,  «   t''*  &Tm«T«r*« 

1.27  C^  *    ' 


1.28       £3  ^ 


T 


where  7/w«  is  proportional  to  the  maximum  vibrational  frequency  1£  , 

c<  is  related  to  the  nearest  neighbor  forces,  and  K  to  the  forces  which 

oppose  the  bending  of  the  chain.   The  approximations  are, 

1.29       f?A)2r  ,      p/K  >  10  }       t~  « 

where  6    is  related  to  the  force  constant  along  the  chain,  and  i      is 
related  to  the  second  nearest  neighbor  forces.   Using  these  relations 
temperature!  regions  are  predicted  for  which  the  heat  capacity  is  at 
first  (lowest  temperatures)  proportional  to  the  cube  of  the  tempera- 
ture, then  changes  to  a  T  dependence,  where  b  lies  between  2.5  and 
3,  than  b  decreases  to  approximately  .8  for  T  about  one  half  Tm. 

To  complete  this  brief  treatment  of  specific  heat  models  which 
might  be  applicable  to  polyethylene,  a  related  phenomenon,  which 
appears  in  the  low  temperature  heat  capacity  of  some  amorphous  solids, 
must  be  mentioned.   If  the  heat  capacity  of  these  substances  is 
measured  calorimetrically  at  low  temperatures,  an  "excess"  specific 
heat  is  observed.   To  understand  the  term  excess  one  should  refer  to 
equation  1.7  which  can  be  rewritten  as  [24 J, 

1.30      Cy  *    5-     13  V       U-l 

where  V  is  the  volume  of  the  sample,  and  W   is  the  velocity  of  sound 
at  these  temperatures,  which  can  be  interpeted  as, 
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1.31     £?  *  T  k  <^<3    ^ 

Therefore  sound  velocity  measurements  should  yield  the  same  value  for 
C/T^  as  those  obtained  from  direct  calorimetric  methods.   However,  it 
has  been  found  experimentally  that  for  some  amorphous  substances,  such 
as  vitreous  germanium,  vitreous  silica,  and  glycerol  glass,  the  calor- 
imetric results  are  much  greater  than  those  given  by  the  acoustic 
methods  and  are  not  T^  dependent  [4,12,17  J. 

The  explanation  for  this  anomolous  effect  is  that  at  low  temper- 
atures there  may  be  modes  excited  which  are  not  measured  by  the 
acoustic  methods.   Since  the  mechanical  sound  measurements  are  con- 
cerned with  the  averaged  (in  the  sense  that  the  measurements  must  be 
macroscopic)  elastic  forces  between  the  particles  making  up  the  solid, 
if  there  are  low  frequency  modes  which  do  not  enter  into  these  elastic 
vibrations,  then  the  calorimetric  measurements  should  be  higher  than 
predicted  by  equation  1.30.   To  explain  the  existence  of  these  extra 
modes,  Rosenstock  has  suggested  that  in  disordered  lattices,  such  as 
the  amorphous  solids,  microscopic  holes  exist  which  may  contain  a 
particle  which  is  bound  in  only  one  or  two  directions  [39,40  J.   These 
loosely  bound  mass  units  can  be  likened  to  particles  free  to  vibrate 
within  a  definite  volume  (the  hole),  and  could  thereby  account  for  the 
low  frequency  non-acoustic  modes.   However,  there  are  other  hypotheses 
which  attempt  to  explain  the  modes  as  being  due  to  other  phenomena 
such  as  non-symmetric  lattice  sites  or  distorted  bonds  []3 ,  27  J .   In 
any  case,  the  origin  of  these  low  frequency  modes  is  not  clear  at 
this  time.   In  the  case  of  the  examples  cited  above  it  has  been  possi- 
ble to  empirically  describe  the  low  temperature  departure  of  the 
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specific   heat    from  the  T3  relation  by  adding  a  term   (or  terms)    pro- 
portional  to   the   Einstein   specific   heat   to   equation   1.30  [2,4,17]    . 
This   implicitly   implys   that    some  of   the  vibrating   units   do   not    con- 
tribute to  the  acoustic  Inodes   and   are  accounted   for    in  the  total 
specific   heat  by  an  expression   such  as, 

1.32  C*x   ~-  2Nt-AB  l&) 

where        -0-£  =    --£■     }     t  ( -qr )  -     \t^h^\)  ) 

V^  is  the  excitation  frequency  for  these  modes,  and  Ms    is  the  num- 
ber of  vibrators  involved  [24 J.   Adding  a  term  such  as  this  to  equation 
1.7  or  1.30  results  in 

rJ  +  ef  (3-) 


1.33 


C  «  A 


where  f\   and  \$   are  constants.   Since  the  specific  heat  is  no  longer 
proportional  to  the  cube  of  the  temperature,  a  plot  of  the  C/T  versus 
temperature  will  demonstrate  the  presence  of  the  excess  heat  capacity. 
2.   EXPERIMENTAL  DETAILS 

Introduction 

In  all  calorimetric  investigations  it  is  necessary  to  know  two 
quantities;  the  temperature  of  the  sample  and  the  acount  of  heat  added 
to  it.   In  this  experiment,  meeting  the  former  requirement  was  the 
more  difficult. 

The  temperature  determination  must  be  made  quickly  and  accurately. 
Since  resistance  thermometry  fulfills  these  requirements,  doped  ger- 
manium resistors  were  selected  as  thermometers.   These  were  chosen 
rather  than  carbon  resistors  since  germanium  resistors  are  reputed 
to  maintain  their  resistance-temperature  relationship  even  after  a 
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large  number  of  cyclings  between  room  temperatures  and  liquid  helium 
temperature  [18,25]  .   Other  experimenters  have  found  that  this  is  not 
true  for  carbon  resistors,  so  that  they  must  be  recalibrated  after 
each  temperature  cycling  [l8]  .   The  need  for  recalibration  would 
greatly  complicate  the  experiment  thus  making  the  choice  of  germanium 
resistors  an  easy  one. 

Since  the  temperature  range  covered  by  this  investigation  was 
from  about  2  to  30°K,  doped  germanium  resistors  (Minneapolis-Honey- 
well, models  2401  and  3401)  were  selected  which  had  sufficiently  high 
sensitivity,  and  yet,  for  the  above  temperatures,  had  resistance 
values  within  the  range  of  ordinary  laboratory  resistance  measuring 
devices.   Once  the  resistors  were  chosen  the  most  important  task  was 
to  find  the  relation  between  the  resistance  and  the  absolute  tempera- 
ture.  The  discussion  of  the  calibration  procedures  is  the  subject  of 
Appendix  A. 

Having  insured  that  a  quick  and  accurate  determination  of  the 
sample  temperature  and  its  changes  could  be  made,  attention  was  di- 
rected to  the  second  requirement.   Adding  heat  electrically  to  a 
resistance  element  in  thermal  contact  with  the  sample  provides  a 
simple  and  accurate  method  of  deliberately  raising  the  sample's  tem- 
perature.  However  great  care  must  be  taken  to  insurfr  that  this 
energy  goes  solely  to  the  sample,  and  that  it  is  the  only  source  of 
heat  input.   To  this  end,  an  adiabatic  calorimeter  was  constructed  in 
which  the  sample  could  be  suspended  in  vaccuo  by  high  thermal  resis- 
tance mountings  from  a  heat  shield,  which  completely  surrounded  the 
sample  and  was  held  at  the  same  temperature.   The  design  and  construc- 
tion of  this  calorimeter  will  be  discussed  in  a  subsequent  section. 
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Sample  Preparation 

The  heat  capacity  of  three  polyethylene  samples  of  various  den- 
sities were  to  be  measured.   Two  of  the  samples  were  Mar lex,  manu- 
factured by  Phillips  Petrol ium  Company,  and  were  received  in  the  form 
of  extruded  cylindrical  pellets  approximately  3  mm.  in  diameter  and 
length.   One  of  the  samples,  Marlex  II,  had  branching  introduced  to 
impede  crystallization.   The  third  sample,  a  typical  commercial  low 
density  polyethylene,  was  a  1"  round  bar.   The  samples  are  referred 
to  as  Marlex  I,  Marlex  II,  and  Low  Density  in  the  remainder  of  this 
thesis. 

An  important  consideration  in  the  preparation  of  the  samples  was 
that  the  time  interval  between  the  application  of  heat  and  the  achieve- 
ment of  thermal  equilibrium  be  as  short  as  possible  so  as  to  minimize 
errors  due  to  heat  leaks  into  or  out  of  the  sample.   The  low  thermal 
conductivity  of  polyethylene  at  low  temperatures,  makes  this  a  diffi- 
cult problem  [ 34].   Because  of  the  calorimeter  design,  the  samples 
were  to  be  cast  into  9/16"  rods  about  6  cm.  in  length.   If  such  a 
polyethylene  rod  were  heated  from  the  center,  it  was  estimated  that 
it  would  require  times  of  about  100  seconds  to  achieve  thermal  equili- 
brium at  5°K  and  nearly  1  1/2  hours  at  30°K.   Such  comparatively  long 
thermal  relaxation  times  would  undoubtedly  introduce  serious  errors 
in  the  determination  of  the  amount  of  heat  added  to  the  sample  unless 
exceptional  precautions  in  the  control  of  the  temperature  of  the  adia- 
batic  shield  were  taken,  in  addition  to  making  the  experiment  extreme- 
ly tedious  to  perform.   Thus  considerable  effort  was  expended  to  pre- 
pare samples  with  shorter  thermal  equilibrium  times.   The  method 
chosen  was  in  inbed  within  the  sample  99.999%  pure  copper  wires  such 
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that  the  wires  were  no  farther  than  .3  cm.  apart.   In  addition  the 
sample  was  wrapped  in  a  copper  foil.   The  heat  could  thus  be  distri- 
buted throughout  the  sample  by  means  of  low  thermal  resistance  copper 
paths.  This  was  expected  to  greatly  reduce  the  thermal  relaxation 
times  estimated  above.   The  success  of  this  method  is  seen  by  the  fact 
that  the  time  required  during  the  actual  experiment  to  achieve  thermal 
equilibrium  was  never  more  than  about  12  minutes  at  the  highest  temp- 
eratures. 

Introduction  of  the  wires  into  the  sample  in  a  uniform  fashion 
presented  a  considerable  problem.   The  following  system  was  eventually 
evolved.   The  pellets  (or  in  the  case  of  the  low  density  sample,  the 
shavings  obtained  from  turning  the  bar  on  a  lathe)  were  first  cast 
into  small  rods  of  1/4"  and  3/16",  and  wrapped  with  30"  of  the  copper 
wire.   The  rods  were  compactly  assembled  within  a  9/16"  mold,  and 
additional  wires  inserted  into  the  resulting  gaps.   The  sample  was 
then  cast  at  155°C  under  900  PSI  for  16  hours,  and  then  allowed  to 
cool  to  room  temperature  (about  12  hours). 

The  extra  material  added  to  the  samples  was  to  be  considered 
part  of  the  calorimeter's  sample  holder.   It  consisted  of  30'  of 
.0051"  diameter  drawn  copper  wire,  .001"  thick  copper  foil  5.580  by 
5.035  cm.,  and  .08  grams  of  G.E.  7031  varnish. 

The  polyethylene  density  was  determined  by  measuring  the  density 
of  the  entire  sample  and  correcting  for  the  copper  wire  which  was 
introduced.   The  crystallinity  of  each  sample  was  then  calculated  by 
use  of  equation  1.1  with  the  values 

2.1     (>A  =  .8-^  cyn/c+»* 


2.2     PQ     '      .  °l   ^  °)        yv*  A**i3 
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where  the  fa  and  Pc   refer  to  pure  amorphous  and  crystalline  poly- 
ethylene densities  Q 9 , 36 J.   It  must  be  noted  here  that  the  densities 
given  above  are  not  the  only  ones  available.   For  the  amorphous  case 
there  is  fair  agreement  in  the  literature  that  .85  is  a  good  average, 
but  more  attention  must  be  given  to  the  crystalline  density  [ 15 J .   For 
this  case,  the  density  has  been  measured  by  X-ray  methods  to  be  as 
high  as  1.014  gm/cm3  and  as  low  as  .964  gm/cm3  depending  on  the  amount 
of  branching  present  [  49  J.   Furthermore  voids  or  free  volume  in  the 
spherulitic  crystalline  form  of  bulk  polyethylene  could  reduce  the 
higher  figure  noted  above  even  in  the  absence  of  branching  [ 10,11  J. 
These  considerations  will  be  referred  to  again  when  the  error  of  the 
experiment  is  discussed.   The  measured  densities  and  calculated  crys- 
tallinities  are  given  in  the  table  below,  where  the  error  indicated 
in  the  latter  quantity  is  that  due  to  the  uncertainty  in  the  density 
measurement  only. 

TABLE  I 
Densities  and  crystallinities  of  the  samples. 

Density  %  Crystal linity 

Marlex  I         .97  3  t    .002  gm/cm3  84.1  t   1.5% 

Marlex  II        .958  t    .002  74.5  t   1.6 

Low  Density      .915  t    .002  45.9  ±   2.5 

Calorimeter  Design  and  Construction 
In  order  to  accurately  measure  heat  capacities,  the  sample  must 
be  thermally  isolated  to  insure  that  the  noted  temperature  differences 
of  the  sample  are  due  entirely  to  the  measured  heat  input.   To  accom- 
plish this,  it  is  necessary  to  remove  all  conduction  and  radiation 
paths  to  and  from  the  sample.   The  method  employed  here  was  to  enclose 
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the  sample  in  a  shield  which  is  kept  at  the  same  temperature  as  the 
sample,  but  insulated  from  it. 

In  this  experiment  the  heat  capacity  would  be  measured  at  various 
temperatures  above  the  temperature  of  a  liquid  helium  bath  by  success- 
ive applications  of  heat  to  the  sample.   Hence  the  shield  could  not  be 
in  good  thermal  contact  with  the  bath.   This  thermal  isolation  was 
accomplished  by  suspending  the  shield  inside  of  a  vacuum  container 
which  was  in  turn  in  contact  with  the  bath.   A  can  type  adiabatic 
shield  was  made  from  a  copper  2"  tube,  3"  long,  threaded  at  the  top 
to  mate  with  its  lid.  The  lid  was  soldered  to  a  3/8"  stainless  tube 
(.01"  wall  thickness),  which  connected  to  a  suitable  vacuum  container. 
The  tube  then  extended  out  of  the  cryostat  to  peripheral  vacuum  equip- 
ment and  served  as  the  pumping  line  and  electrical  lead  housing.   The 
shield  also  had  provisions  for  receiving  a  resistor  which  would  sense 
the  shield's  temperature. 

There  would  be  a  heat  leak  from  the  shield  to  the  vacuum  con- 
tainer through  the  support,  but  this  was  compensated  for  by  means  of 
an  electrical  heater  on  the  shield.   In  addition  this  heater  was  used 
to  raise  the  shield's  temperature  to  that  of  the  sample.   Seventeen 
feet  of  .002"  manganin  wire,  coiled  and  glued  with  G.E.  7031  varnish 
to  the  top  of  the  shield,  served  as  the  shield  heater. 

The  sample  holder,  in  addition  to  holding  the  9/16"  sample,  was 
used  to  mount  the  sample  heater,  temperature  sensing  resistor,  and  a 
receptacle  which  would  be  the  contact  point  for  the  thermal  switch 
(to  be  discussed  later).   A  copper  ring  with  the  necessary  provisions 
was  constructed  for  this  purpose.   When  heat  capacity  measurements 
were  made,  the  measured  values  would  be  the  total  of  the  sample  plus 
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the  addenda  (the  holder,  heater,  thermometer,  etc.).   Although  the 
latter' s  heat  capacity  would  be  determined  experimentally  and  sub- 
tracted from  the  total,  it  was  still  desirable  to  keep  it  as  small 
as  possible.   Based  on  the  size  of  the  calorimeter,  it  was  determined 
that  a  sample  of  10  cm^  could  be  accommodated.   With  a  sample  of  this 
size  it  was  estimated  (correctly,  as  was  determined  later)  that  the 
heat  capacity  of  the  addenda  would  amount  to  less  than  10%  of  the 
total,  which  was  deemed  acceptable. 

The  sample  heater,  which  consisted  of  15"  of  .002"  manganin  wire 
(100  ohms),  was  wouftd  around  the  thermometer  receptacle.   A  four  wire 
connection  was  made  to  the  heater,  and  the  wires  were  thermally 
grounded  to  the  adiabatic  shield.   The  power  was  measured  potent io- 
metrically  using  a  Leeds  and  Northrup  K-3  potentiometer.   A  Varian 
model  501A  time  mark  generator,  used  in  conjunction  with  a  Hewlett- 
Packard  model  521C  electronic  counter,  determined  the  heating  times 
(the  timer  was  later  changed  to  an  Anadex  model  CF-200R  counter- 
timer)  and  was  triggered  by  the  sample  heater  switch. 

It  would  be  impossible,  regardless  of  the  care  taken,  to  insure 
that  the  shield  would  be  at  exactly  the  same  temperature  as  the  sample 
at  all  times.   However,  if  the  shield  were  thermally  isolated  from 
the  sample,  small  temperature  differences  could  be  tolerated.   Based 
on  this  the  sample  holder  was  suspended  within  the  shield  by  6  nylon 
threads,  approximately  .12  mm.  in  diameter,  and  a  minimum  of  1.9  cm. 
in  length.   Assuming  a  maximum  of  .5°  temperature  difference  between 
the  shield  and  the  sample,  the  heat  leak  through  these  supports  would 
be  on  the  order  of  10"-3  mj'/sec  which,  for  equilibrium  times  of  less 
than  an  hour,  could  be  neglected  (radiation  losses  are  of  the  same 
order  of  magnitude)  . 

30 


When  the  calorimeter  was  initially  cooled,  or  when  cooling  back 
to  the  bath  temperature  after  a  series  of  measurements,  a  method  had 
to  be  incorporated  in  the  design  to  cool  the  isolated  sample  in  a 
reasonable  length  of  time.   For  this  reason  a  mechanical  switching 
device  was  constructed.   This  was  chosen  over  the  use  of  exchange  gas 
in  the  vacuum  spaces,  since  the  latter  method  has  been  shown  to  have 
quite  long  recovery  times  and  has  attendent  gas  absorption  problems 
[29j.   The  switch  consisted  of  a  brass  finger  which  would  make  contact 
with  a  receptacle  on  the  sample  holder.   The  finger  was  soldered  to  a 
brass  rod  which  extended  out  of  the  vacuum  container  into  the  liquid 
helium  bath  through  a  nickel  bellows  and  from  there  out  of  the  cryo- 
stat.  The  finger  could  be  raised  and  lowered  manually  so  as  to  make 
thermal  and  mechanical  contact  with  the  holder  whenever  it  was  re- 
quired.  The  original  design  was  modified  after  an  initial  trial  by 
tinning  the  brass  finger  and  its  mating  receptacle  with  metallic  in- 
dium, for  better  thermal  contact.   In  addition  a  brass  leaf  was  added 
which  would  make  contact  with  the  adiabatic  shield  when  the  thermal 
switch  was  closed  in  order  to  cool  the  entire  calorimeter  simultane- 
ously. 

The  two  thermometers  (resistors)  were  attached  to  the  shield  and 
the  sample.   Three  wire  connections,  which  were  thermally  anchored  to 
the  shield,  were  made  to  each  resistor.   The  resistance  measuring 
equipment  was  identical  to  that  used  by  Reese  and  Tucker  and  will  not 
be  discussed  here  [34 J. 

The  completed  calorimeter  is  shown  in  figure  1. 
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Figure  1,      Calorimeter  with   sample   in  place 


a.  "*  Triple 

b.  Samole  holder 

c.  Sample  thermometer 

d.  Sample  heater 

e.  Thermal  switch' 

f.  Shield  thermometer 
■  .  Nylon  support 

h .  Anchor 


Measurement  Procedures 

In  preparation  for  the  heat  capacity  measurements,  the  sample 
was  placed  in  the  calorimeter  and  cooled  to  the  desired  initial  tem- 
perature with  the  thermal  switch  closed.   When  the  thermometers  indi- 
cated that  the  sample  and  shield  were  at  the  liquid  helium  bath  temp- 
erature, the  thermal  switch  was  then  opened.   This  set  up  a  small 
amount  of  frictional  heating  as  the  contact  was  broken,  but  within 
30  minutes  the  system  would  again  be  in  equilibrium. 

The  sample  was  then  heated  by  applying  electrical  power  to  the 
sample  heater  in  amounts  which  would  raise  the  temperature  of  the 
sample  in  5%  increments  of  its  average  temperature.   The  power  was 
recorded  at  each  temperature  step  as  well  as  the  initial  and  final 
readings  of  the  time  and  the  thermometers.   The  heating  times  were 
kept  between  20  and  60  seconds,  which  required  that  the  power  be  in- 
creased every  four  to  six  measurements.   When  the  desired  final  temp- 
erature was  reached,  the  thermal  switch  was  again  closed,  the  sample 
cooled,  and  the  procedure  repeated. 

While  the  heat  was  being  supplied  to  the  sample,  the  adiabatic 
shield  was  also  heated  by  means  of  its  heater.   This  required  care 
since,  as  previously  mentioned,  it  was  highly  desirable  that  the 
shield  and  the  sample  be  at  the  same  temperature.   After  several 
trials,  however,  certain  regularities  in  the  amount  of  heat  needed 
for  the  shield  at  various  temperatures  became  evident,  and  the  desired 
result  was  obtained.   It  is  felt  that  there  never  existed  more  than  a 
one  percent  temperature  difference  between  the  two  thermometers  while 
actual  measurements  were,  being  taken,  and  for  the  bulk  of  the  runs, 
the  difference  was  considerably  less. 
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Due  to  the  close  proximity  of  sample  thermometer  to  the  heater, 
while  the  power  was  being  supplied  the  temperature  would  be  observed 
to  overshoot  the  desired  value,  and,  when  the  power  was  removed,  drop 
back  exponentially.   The  time  required  for  the  sample  to  reach  equil- 
ibrium was  from  2  seconds  to  about  12  minutes  as  the  temperature  was 
increased  from  its  lowest  value  to  its  highest.   These  times  were  not 
considered  excessive,  except  at  the  highest  temperatures  when  the 
shield  might  tend  to  drift  while  equilibrium  was  being  established. 
This  drifting  was  countered  by  monitoring  both  thermometers  and  keep- 
ing the  shield  at  the  desired  temperature  by  removing  or  increasing 
the  shield  power. 

The  heat  capacity  of  the  samples  was  determined  by  five  separate 
sets  of  measurements.   The  first  three  covered  the  range  of  4.5  to 
30°K,  one  series  for  each  of  the  three  samples.   The  other  two  were 
for  the  temperature  interval  2.4  to  about  8°K  for  the  Marlex  I  and 
Low  Density  samples  only.   For  each  set  of  measurements  the  aforemen- 
tioned temperature  ranges  were  covered  three  times. 

Since  the  measured  heat  capacities  were  the  sum  of  those  of  the 
sample  and  the  addenda,  the  latter  was  measured  in  two  separate  ex- 
periments, for  the  temperature  ranges  of  4.5  to  30°K  and  2.4  to  8°K, 
so  that  its  contribution  could  be  subtracted  from  the  total.   In  each 
of  the  experiments,  the  addenda  heat  capacity  was  measured  four  times 
over  the  covered  range.   The  procedures  were  identical  to  those  des- 
cribed above  except  due  to  the  reduced  mass,  the  amount  of  heat  re- 
quired and  the  equilibrium  times  were  both  substantially  decreased. 
3.   RESULTS  OF  THE  EXPERIMENT 

Data  Reduction 

When  the  experimental  measurements  of  the  heat  capacity  were 
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made,  the  data  consisted  of  the  power  reading,  the  elapsed  time  inter- 
val during  which  the  power  was  applied,  and  the  initial  and  final 
temperatures  as  indicated  by  resistance  values.   The  heat  capacity- 
can  then  be  calculated  from, 

3.1  C  =   Power  X  Elapsed  Time . 

TFinal  "  TInitial 

Due  to  the  complexity  of  the  resistance-temperature  relation  (Appen- 
dix A),  the  calculations  were  performed  by  the  CDC  1604  digital  com- 
puter. 

All  of  the  data  were  reduced  to  heat  capacities,  and  of  the  ones 
which  displayed  large  deviations,  only  those  attributable  to  obvious 
blunders  in  data  recording  were  rejected. 

Since  the  heat  capacities  measurements  for  each  sample  represen- 
ted the  total  of  the  sample  plus  the  addenda,  the  latter' s  contribu- 
tion had  to  be  taken  into  account  and  corrected  for.   As  mentioned 
before,  measurements  of  the  heat  capacity  of  the  addenda  were  per- 
formed.  These  heat  capacities  were  for  discrete,  random  temperatures 
and,  as  such,  were  not  useful  for  the  adjustment  which  was  to  be  made. 
Therefore  the  data  were  best  fitted,  in  the  least  squared  sense,  by 
the  computer  to  a  continuous  heat  capacity- temperature  relation  of 
the  form, 

3.2  C  =  ^i  ^T 

With  this  relationship  the  total  heat  capacity  for  eacli  sample 
could  be  corrected  so  that  the  specific  heat  of  the  sample  alone 
could  be  determined.   These  calculations  were  also  performed  by  the 
computer  and  the  resulting  heat  capacity  data  were  again  fitted  to 
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equation  3.2.   The  specific  heat  was  obtained  by  dividing  the  heat 
capacity  by  the  mass  of  the  sample. 

The  measured  specific  heat  for  each  sample  is  found  in  Appendix 
C,  and  the  parameters  of  equation  3.2  for  the  samples  and  the  addenda 
in  Appendix  E.   The  smoothed  results  from  the  fitted  curve  are  shown 
in  Appendix  D,  and  a  specific  heat -temperature  graph  is  displayed  in 
figure  2  where  the  difference  between  the  smoothed  curve  and  the  actual 
data  points  is  too  small  to  be  shown. 

The  temperature  range,  total  number  of  data  points,  and  the  per 
cent  RMS  deviation  of  the  actual  data  from  the  smooth  curve  are  pre- 
sented in  Table  II. 

Table  II 

Temperature  range,  amount  of  data,  and  RMS  fit  of  the  data 
to  a  smoothed  curve  for  each  of  the  samples  and  the  addenda. 


Temperature 
range 

Number  of 
data  points 

%RMS 
deviation 

Addenda 

2.4  -    30°K 

218 

1.70% 

Marlex  I 

2.4  -    30 

161 

1.52 

Marlex  II 

4.5   -    30 

117 

1.01 

Low  Density 

2.4  -    30 

158 

1.09 

Extrapolation  to  100%  and  0%  Crystallinity 
As  mentioned  in  the  first  chapter,  the  specific  heat  of  any 
polyethylene  sample  might  best  be  expressed  as  admixture  of  the  specif- 
ic heats  of  the  completely  amorphous  and  completely  crystalline  phases. 
As  suggested  by  the  earlier  low  precision  work  of  Reese  and  Tucker,  and 
very  well  substantiated  by  the  present  investigation,  the  relation  is 
linear  [34].   That  is 

3.3     C  =  XCI00    +(l-*)C0 
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Figure  2.   Specific  heat  versus  temperature. 
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where  the  C's  on  the  right  side  refer  to  hypothetical  values  for  the 
specific  heat  of  crystalline  and  amorphous  polyethylene.   This  lin- 
earity is  reasonable  if  one  considers  that  it  is  possible  (in  theory) 
to  break  a  sample  up  into  many  very  small  macroscopic  volumes,  each 
of  which  is  either  entirely  amorphous  or  crystalline  (if  the  boundar- 
ies are  neglected).   The  two  idealized  polyethylenes  should  have  dif- 
ferent values  for  the  specific  heat  at  any  particular  low  temperature* 
since  the  amorphous  state  should  be  characterized  by  a  set  of  para- 
meters which  will  differ  from  those  of  the  crystalline  state  inasmuch 
as  the  interchain  binding  scheme  is  completely  different. 

Using  this  assumption,  the  heat  capacities  for  the  three  sam- 
ples, divided  by  the  cube  of  the  temperature,  were  plotted  as  a 
function  of  crystallinity  for  various  temperatures  (Figure  3).   The 
linearity  is  quite  apparent,  and  on  the  basis  of  this,  it  is  felt  that 
equation  3.3  is  indeed  valid.   Using  this  linear  relationship  one  can 
perform  a  simple  extrapolation  to  obtain  values  for  the  specific  heat 
as  a  function  of  temperature  for  the  hypothetical  1007o  and  0%  cases. 
However,  rather  than  use  the  C/T3  relation,  which  is  shown  for  des- 
criptive purposes  only,  the  computer  was  used  to  perform  the  extrapo- 
lation from  the  equation 

3.4     C  =  W  *  ~t"  -v" 

where  /vv\  and  Xr  were  determined  from  the  experimental  (smoothed) 
values  of  C,  the  specific  heat  for  the  three  samples  of  differing 
crystallinities .  The  extrapolation  to  both  0%  and  100%  was  performed 
every  .04°  and  the  results  were  fitted  to  equation  3.2,  the  same  as 
the  real  samples. 

Appendix  D  presents  the  extrapolated  data,  and  figure  2  is  a 
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plot  of  the  specific  heat  versus  temperature.   In  Appendix  E  are 
shown  the  parameters  for  equation  3.2  for  all  of  the  real  and  extra- 
polated samples. 

The  100%  Extrapolated  Data 
The  heat  capacity  obtained  from  the  extrapolation  to  100% 
crystal linity  had  some  interesting  features.   A  plot  of  C/T^  versus 
temperature  shows  a  relatively  large  temperature  region  (up  to  about 
9°K)  where  the  specific  heat  has  a  cubic  temperature  dependence  (fig- 
ure 4).   This  is  in  agreement  with  the  expectations  if  there  exist 
relatively  large  binding  forces  in  the  crystalline  case,  since  on  the 
basis  of  the  Debye  model  larger  forces  imply  larger  & S  which  in  turn 
imply  an  increase  in  the  temperature  range  in  which  the  cubic  depen- 
dence holds  (equation  1.7). 

For  this  case,  an  excess  heat  capacity  is  not  indicated  by  the 

3 
C/T  plot,  nor  is  one  expected  since  this  material  should  be  a  con- 
tinuous polycrystal. 

When  the  heat  capacity  for  this  hypothetical  sample  is  plotted 
against  temperature  on  log-log  paper,  it  is  possible  to  obtain  the 
parameters  A  and  a  of  the  expression, 

3.5    C=ATa. 

There  appeared  two  regions  in  which  they  were  constant,  indicating 
that  the  specific  heat  is  proportional  to  some  power,  a,  of  the  temp- 
erature. The  first  region  was  the  cubic  dependence  noted  before.   The 
other  was  from  about  20  to  26°K  where  the  exponent,  a,  was  approximate- 
ly 2.37-2.40.   In  between  the  value  decreased  uniformly.   This  is  in 
disagreement  with  the  prediction  of  Genensky  and  Newell  who  expected 
it  to  be  equal  to  some  constant  slightly  greater  than  2.5  for  an 
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Figure  4. 
A  plot  of  C/T3  versus  temperature  for  polyethylene, 
where  the  solid  lines  are  the  experimental  and  the 
extrapolated  data,  the  dotted  line  below  the  0%   curve 
is  the  data  minus  the  Einstein  term,  and  the  other 
dotted  lines  are  the  predicted  T  =  0°K  values  for 
these  samples.  The  dots  represent  the  predicted 
specific  heat  values  in  all  cases. 
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appreciable  temperature  region,  for  a  crystalline  polymer  [l9J. 
This  contradiction  will  be  discussed  later. 

The  Tarasov  model  could  be  well  fitted  to  the  extrapolated  data. 
The  method  used  to  determine  the  parameters  was  essentially  the  same 
as  that  suggested  by  Wunderlich,  except  that  the  equations  used  are 
for  the  extended  model,  equation  1.19  rather  than  the  one  suggested  by 
Tarasov,  which  doesn't  separate  the  longitudinal  and  transverse  con- 
tributions [52].   According  to  equation  1.22  there  should  exist  a 
region  in  which  the  "heat  capacity  becomes  linear  in  temperature.   This 
will  occur  at  a  temperature  well  above  that  at  which  most  of  the  three 
dimensional  modes  are  filled  and  above  the  temperature  range  of  this 
investigation.   Therefore  the  data  of  Dainton,  et  al,  giving  the 
specific  heat  of  polyethylene  from  20-300°K,  were  used  to  determine 
the  linear  region  [l3j.   This  procedure  will  be  justified  later.   The 
linear  region  appears  to  occur  between  80  and  110°K,  where  C/T  has 
the  value  6. 69mj/°K  -gm.   Rewriting  equation  1.22,  it  is  seen  that  in 
this  region, 

3.6      C/T  =  7T1R[^(-^Tt  +"^i.)I- 

Another  relation  between  "®"it  and  "8^1-  is  necessary  to  determine  them 
uniquely.   This  was  accomplished  by  considering  that  the  "Oin  since  it 
was  the  longitudinal  mode  limit  for  the  chain,  might  be  well  character- 
ized by  the  largest  frequency  for  which  the  C-C  bond  will  be  stretched. 
Reference  to  the  Raman  spectrum  of  polyethylene  indicates  a  value  of 
1150  cm"1  for  this  parameter,  and  therefore  a  value  of  1655°K  was 
assumed  for  ■©'it  [l5].   This  fixed  -#It  at  707°K.   (The  use  of  an  appar- 
ent optical  mode  in  this  discussion,  which  up  to  now  had  only  consid- 
ered acoustical  vibrations,  will  be  discussed  later). 
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Having  determined  these  two  parameters,  it  is  then  necessary  to 
fix  the  remaining  two  thetas  needed  for  the  extended  Tarasov  model, 
-£}r  and-^c  •   One  relation  comes  from  equation  1.20,  which  is  here 
rewritten, 

3>7       C/T3   =  ZR  *p[~k  (*«*&'*  3}I3aJ]  ' 

This  is  the  low  temperature  limit,  where  C/T^  is  constant,  and  the  re- 
sults below  r^  9°K  could  be  used  to  evaluate  the  left  hand  side.   In 
order  to  solve  for  either  "^ir  or  ~&3L    another  relation  is  needed. 
This  was  obtained  from  sound  velocity  measurements  at  liquid  nitro- 
gen temperatures  for  a  highly  crystalline  sample.   These  measurements 
were  performed  by  Professor  W.  Reese  [ 33  J.   Extrapolation  of  the  re- 
sults to  zero  degrees  Kelvin  predicted  that  the  ratio  of  the  longi- 
tudinal sound  velocity  to  that  of  the  transverse  was  1.67.   Since  the 
thetas  have  the  same  relative  relationship  to  one  another  as  do  the 
sound  velocities,  "^Jr  and  -&31.   could  now  be  determined.  The  Tarasov 
parameters  for  this  case  are  found  in  Table  III. 

The  0%  Extrapolated  Data 
Extrapolation  to  the  0%  crystallinity  (completely  amorphous) 
case  led  to  results  which  were  quite  different  from  those  discussed 
above.   In  particular,  the  plot  of  C/T3  versus  the  temperature  shows 
that  the  heat  capacity  was  not  proportional  to  the  cube  of  the  temp- 
erature anywhere.   The  curve  starts  with  a  positive  slope  at  2.5°K, 
increases  to  a  maximum  at  around  5°K,  then  decreases  uniformly.   This 
is  also  the  case  for  the  real  samples  as  can  be  seen  in  figure  4. 
Since  even  the  Marlex  I  sample  (crystallinity  of  84.1%)  displays  this 
hump,  while  the  1007o  extrapolated  sample  does  not,  it  is  quite  apparent 
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that  the  presence  of  the  hump  must  be  due  to  the  amorphous  poly- 
ethylene.  The  hump  is  interpeted  as  an  indication  of  an  excess  heat 
capacity.   Referring  to  the  first  section  the  heat  capacity  for  the 
lowest  temperatures  should  then  be  of  the  form, 


C  =AT^  -V^f^ 


3.8 

,.,        y-  3J.M 

where  <>,  ^  and  ^  are  constants  to  be  determined.   This  equation 
was  best  fitted,  in  the  least  squared  sense,  to  the  extrapolated  re- 
sults (and  to  the  Low  Density  data  where  a  correction  was  applied  to 
account  for  the  crystalline  part)  to  yield  values  for  o  ,   A     and  ~&£  . 
For  the  model  which  describes  the  portion  of  the  heat  capacity  not 
given  by  the  Einstein  term,  one  is  forced  to  accept  a  continuum  type, 
since  it  is  not  possible  to  do  lattice  dynamics  for  the  amorphous 
case.   From  the  foregoing  discussion,  the  extended  Tarasov  model  was 
chosen.   Assuming  this,  A   is  interpeted  as 


3.10 


A  ;    dJL*(3fi  -  *)  ["3  (  *,r  35  +  *>..*£  )J  • 


There  is  no  reason  to  assume  that  the  one  dimension  parameters 
will  change  from  those  used  previously,  since  they  are  merely  func- 
tions of  the  forces  along  the  chain.   However  another  relation  be- 
tween wjf  and  ^t  is  still  needed  in  order  to  determine  them  uniquely. 
The  ratio  of  the  longitudinal  and  transverse  sound  velocities  used 
earlier  is  not  really  applicable  in  this  case,  since  the  sound  vel- 
ocities were  those  for  a  highly  crystalline  sample;  but,  lacking  any 
other  evidence,  the  same  relation  was  assumed. 

Using  the  above  assumptions,  and  the  values  of  r\   and  0  ,  the 
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Tarasov  parameters  can  be  found.   They  are  presented  in  Table  III 
along  with  other  results  for  this  extrapolated  hypothetical  sample. 

Summary 

The  percentage  difference  between  the  0%  and  the  100%  extrapola- 
ted samples  were  plotted  against  temperature  as  shown  in  figure  5. 
Although  the  curve  is  approaching  zero,  it  is  impossible  to  determine 
with  any  degree  of  certainty  when  if  ever  this  would  occur.   However, 
it  is  felt  that  the  percentage  difference  would  be  insignificant  by 
100°K,  where  the  data  is  linear  in  temperature  [l3j.   This  supports 
the  assumption  that  "£-,  will  be  the  same  for  all  samples. 

The  heat  capacities  for  the  real  samples  and  the  extrapolated 
ones  are  shown  in  Appendix  D  at  .  5°K  increments,  together  with  the 
specific  heats  predicted  by  the  Tarasov  equation,  equation  1.19,  which 
is  here  expanded  as 


3.11 


+r,.x)(,_X)[a(D,^)-^[p,(^-D3(^)]] 


The  superscripts  on  the  ~&S  refer  to  either  the  0%  or  the  100%  para- 
meters,  which  are  presented  in  Table  III.   The  crystallinities,  /   , 
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are  shown  in  Table  I. 

Two  hitherto  undefined  terms  also  appear  in  Table  III.   These 
are  the  Tarasov  parameters  which  would  result  if  no  separation  into 
the  longitudinal  and  transverse  modes  were  assumed,  which  is  the  ori- 
ginal model  discussed  by  Tarasov.   In  this  case,  equation  3.11  would 
become, 

,.         COR(x[o,(i>^[D,(f)-D^l] 


In  addition,  the  classic  Debye  temperature  is  shown  as  calculated 
from  the  low  temperature  limit. 

In  figure  4  the  dotted  line  beneath  the  0%  crystallinity  sample 
curve  is  the  result  of  subtracting  the  term  tL{y   J  f  rom  the  data,  to 
illustrate  how  it  accounts  for  the  excess  heat  capacity.   For  the 
other  samples,  the  dotted  lines  indicate  the  0  K   intercept,  and  for 
all  cases,  the  circled  points  are  the  values  predicted  from  equation 
3.11. 

In  Appendix  D  the  smoothed  values  of  the  specific  heat  for  each 
sample  is  shown,  along  with  the  calculated  specific  heat  of  the  Tara- 
sov model  as  given  by  equation  3.11.   This  calculation  was  performed 
by  the  computer  (the  mass  of  the  polyethylene  vibrating  unit  is  taken 
as  14.1  atomic  mass  units). 
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Table  III 

The  Debye,  Tarasov,  and  excess  heat  capacity  parameters, 
and  the  low  temperature  values  of  the  specific  heat  for 
the  extrapolated  data. 


100% 

0% 

t 

(mj/° 

K-gm) 

0.0 

3.24 

<% 

(°K) 

260.1 

175.4 

-&, 

(°K) 

873.9 

873.9 

*<1 

(°K) 

141.9 

78.6 

-GL 

(°K) 

1655.0 

1655.0 

■Oh 

(°K) 

707.0 

707.0 

■0-3L 

(°K) 

222.8 

123.4 

^3T 

(°K) 

133.7 

74.1 

•6* 

(°K) 

- 

23.0 

c-*a(%) 

(mi/°] 

K-gm,    low 

.00787 

.0256 

"J"*3    temp.) 

Accuracy  of  the  Results 

To  estimate  the  accuracy  of  the  heat  capacity  data,  four  sources 
of  error  were  considered;  the  determination  of  the  amount  of  heat  in- 
troduced into  the  sample,  the  temperature  determination,  the  effect 
of  the  shield,  and  the  heat  capacity  of  the  addenda. 

The  amount  of  heat  energy  introduced  into  the  sample  was  obtained 
from  three  measurements;  voltage,  current,  and  time.   The  voltages  (and 
current)  could  be  determined  to  six  places  on  the  potentiometer.   Since 
the  same  power  was  used  on  from  four  to  six  successive  measurements, 
it  was  possible  to  check  the  reproducibility  of  the  readings.   Only 
at  the  very  highest  temperatures,  where  the  resistance  of  the  heater 
was  probably  changing  with  temperature,  did  it  vary  as  much  as  .1%. 
The  current  measurement  had  the  same  inconsistency,  but  since  it  was 
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determined  by  measuring  the  voltage  across  a  standard  resistor,  it 
also  had  an  inherent  maximum  error  of  .05%.   Timing  was  presumed  ac- 
curate to  one  part  in  thirty  thousand.   The  maximum  total  error  from 
this  source  is  therefore  estimated  to  be  .05%  at  the  low  end  of  the 
temperature  sQSle  and  .25%  at  the  high  end. 

As  discussed  in  Appendix  A,  the  thermometers  were  assumed  cor- 
rect to  within  .5%  of  the  absolute  temperature.   This  uncertainty  will 
influence  the  accuracy  of  the  heat  capacity  results  in  two  ways.   Since 
a  temperature  difference  is  divided  into  the  total  heat  to  obtain  the 
heat  capacity,  an  error  can  arise  there.   However  any  experimental  or 
systematic  errors  in  the  temperature  determination  will  be  the  same 
over  a  short  temperature  interval  and  should  cancel  out  (this  is  en- 
sured by  the  use  of  the  smooth  fitted  curve  which  should  average  out 
these  errors).   Therefore  the  possible  error  from  this  source  should 
be  only  that  due  to  the  error  of  judgment  in  reading  the  resistance 
bridge.   This  results  in  uncertainties  in  the  heat  capacity  determi- 
nation of  .12%  at  2.4°K  up  to  about  .2%  at  30°K. 

The  other  manner  in  which  the  inaccuracy  of  the  temperature 
scale  enters  into  the  error  is  in  assigning  a  temperature  to  a  par- 
ticular heat  capacity  measurement.   Since  the  specific  heat  can  be 
expressed  at  any  point  as  being  proportional  to  some  power  of  the 
temperature  (the  exponent  will  be  between  3  and  2.4  for  this  case), 
the  percentage  error  in  the  specific  heat  will  be  larger  than  that  of 
the  temperature.   In  this  case  the  maximum  possible  error  amounts  to 
1.5%  at  the  lowest  temperature  and  1.2%  at  the  highest  from  this 
source. 

Since  it  was  impossible  to  completely  isolate  the  sample  from 


48 


the  shield  thermally,  considerable  care  had  to  be  taken  to  ensure 
that  they  were  at  the  same  temperature,  to  prevent  possible  heat 
leakage.   Experimentally,  the  effects  of  large  temperature  differ- 
ences were  investigated,  by  purposely  overshooting  the  shield  tem- 
perature and  observing  the  effect  it  had  on  heat  capacity  data.   On 
the  basis  of  these  experiments  it  was  decided  that,  with  the  experi- 
mental procedures  used  when  actual  data  were  being  taken,  this  con- 
tribution to  the  error  could  be  neglected  for  all  but  the  highest 
temperatures,  when,  due  to  the  insensitivity  of  the  shield  thermom- 
eter, difficulty  insuring  equal  temperature  of  the  shield  and  sample 
might  occur.   At  these  temperatures,  the  error  might  approach  .5%. 

The  heat  capacity  of  the  addenda  was  measured  over  the  entire 
temperature  range,  so  that  it  could  be  applied  as  a  correction  to 
the  total  measured  for  any  sample.   Since  the  mass  of  the  addenda 
was  chosen  to  be  as  small  as  possible,  so  that  its  heat  capacity 
would  be  on  the  order  of  10%  of  that  of  the  samples,  any  possible 
errors  which  are  present  in  the  measurement  of  its  heat  capacity 
would  be  reduced  by  a  factor  of  ten  when  applied  to  the  total  possi- 
ble error  for  a  sample.   Considering  the  preceeding  this  uncertainty 
amounts  to  approximately  .2%  over  the  entire  temperature  range. 

Considering  all  of  the  above,  it  is  felt  that  the  measured  heat 
capacities  were  accurate  to  within  1.9%  for  the  temperature  up  to 
about  26°K,  and  3%  above  that. 

The  subject  of  Appendix  B  is  the  measurement  of  the  specific 
heat  of  copper.  The  results  are  compared  with  data  found  in  the 
literature  which  presumably  is  of  higher  precision  and  accuracy  than 
this  work.   The  comparison  shows  that  the  results  are  consistant, 
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within  the  error  estimate,  but  the  present  results  tend  to  be  low 
over  most  of  the  temperature  range,  which  could  be  an  indication  of 
a  systematic  error.   This  point  was  not  investigated  further. 

As  mentioned  before,  there  was  no  statistical  criteria  estab- 
lished for  rejecting  any  of  the  data,  and  the  only  points  discarded 
were  ones  in  which  a  large  deviation  could  be  attributed  to  a  blunder 
in  data  accrual. 

When  the  data  were  fitted  to  equation  3.2,  the  RMS  deviation 
(Table  II)  was  found  to  be  well  within  the  estimated  accuracy  of  the 
experiment  and  therefore  the  use  of  the  smoothed  data  is  considered 
to  be  a  valid  representation  of  the  results. 

Since  one  of  the  objectives  of  this  investigation  was  to  verify 
the  density  (crystal linity)  dependence  of  the  specific  heat  of  poly- 
ethylene, and  an  extrapolation  was  made  to  the  hypothetical  100%  and 
0%  crystalline  phases,  an  estimate  was  made  as  to  the  reliability  of 
this.   The  heat  capacity  of  the  three  real  samples  would  be  exactly 
linear  in  crystallinity  over  the  entire  temperature  range,  if  a  varia- 
tion of  at  most  4%  were  made  in  any  one  of  the  heat  capacities,  or  if 
a  1.5%  adjustment  were  made  in  one  of  the  crystal linities.   Addition- 
ally there  are  other  combinations  of  adjustments  which  could  be  made, 
but  it  is  the  uncertainty  in  the  degree  of  crystallinity  which  is  be- 
lieved to  be  subject  to  the  most  error.   As  noted  in  Table  I,  there 
is  an  inherent,  approximate  2%  error  which  is  due  solely  to  the  un- 
certainties in  the  densities.   In  addition,  since  the  densities  of 
the  amorphous  and  crystalline  phases  are  not  well  established,  it  is 
estimated  that  the  crystallinity  of  a  sample  could  be  in  error  as 
much  as  4%, abased  on  the  density  method  of  determining  crystallinities 
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(equation  1.1).   In  view  of  this,  it  is  felt  that  the  linearity  is 
valid  to  the  precision  of  this  experiment. 

Due  to  the  possible  error  in  the  crystal 1 inities ,  the  accuracy 
of  the  extrapolated  samples  can  be  better  than  2%  for  the  amorphous 
case,  and  .8%  for  the  crystalline  phase.   The  uncertainties  which  are 
present  in  the  actual  specific  heat  is  taken  over  into  the  extrapola- 
ted data.   If  the  relative  error  in  the  heat  capacity  at  any  particular 
temperature  is  the  same  for  all  of  the  samples,  as  should  be  expected 
since  they  should  exhibit  the  same  systematic  error  if  it  is  present, 
the  maximum  possible  error  in  the  extrapolated  data  due  to  this  cause 
will  be  the  same,  or  about  1.9%  over  most  of  the  temperature  range. 
Thus  the  uncertainty  in  these  figures  will  be  approximately  3  and  4% 
for  the  completely  crystalline  and  amorphous  phases. 

The  precision  of  this  experiment  can  be  estimated  by  the  close- 
ness of  the  fit  of  the  actual  data  points  to  the  smoothed  curve  given 
by  equation  3.2.   The  deviations  appear  random  and  the  per  cent  RMS 
deviation  is  given  in  Table  II,  where  the  %  RMS  deviation  is  defined 
by, 

3.13      %RMS  =  100%   I  Hi J   . 

This  precision,  which  is  also  an  estimate  of  the  reproducibility 
of  the  results,  is  between  1.01  and  1.52%  for  the  polyethylene  sam- 
ples. 
4.   CONCLUSIONS 

Comparison  with  Other  Results 

The  results  of  this  research  have  been  compared  with  all  of  the 
known  experimental  work  on  the  specific  heat  measurements  on  poly- 
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ethylene  in  their  common  temperature  ranges.   A  detailed  survey  of 
these  comparisons  follows: 

(a)  Isaacs  and  Garland  [  21 ] 

The  specific  heat  measurements  of  Isaacs  and  Garland  were  on 
a  single  high  density  sample  in  the  temperature  range  1.8  to  5.3°K. 
They  assumed  that  the  specific  heat  in  this  range  displaced  T^  depen- 
dence and  therefore  did  not  separately  determine  the  heat  capacity  of 
the  addenda,  but  rather  calculated  it  and  subtracted  it  from  the  total. 
This  led  as  they  report  to  an  unaccountable  deivation  from  the  T  be- 
havior, but  they  attributed  it  to  experimental  error.   However  it  is 
felt  that  they  were  in  reality  observing  the  excess  heat  capacity 
which  was  disguised  by  their  method  of  analysis. 

For  the  comparison  in  the  mutual  temperature  range,  it  was  noted 
that  both  the  amorphous  and  crystalline  parts  of  polyethylene  are  pro- 
portional to  the  cube  of  the  temperature  if  the  Einstein  term  is  sub- 
tracted.  Since  their  sample  did  not  have  the  same  degree  of  crystallin- 
ity  of  any  of  the  ones  measured  in  this  work,  the  following  equation 
was  used  to  determine  what  specific  heat  should  be  expected  on  the 
basis  of  the  present  results; 

This  is  merely  the  low  temperature  limit  of  equation  3.12  with  the 
measured  values  inserted.   This  assumed  no  model,  however;  only  that 
the  specific  heat  be  T^  dependent  at  low  temperatures,  which  it  is 
for  any  realistic  case,  and  that  the  linearity  in  crystallinity  is 
valid.   For  the  crystallinity  of  Isaacs  and  Garland's  sample,  two 
values  were  investigated;  first  the  reported  X-ray  determined  result 


52 


of   90%,    and   second,    a  crystal  Unity  of   85%  deduced   from  the  reported 
density.      In  Table   IV  where  the  comparison   is  made,    the  two   different 
values   for   the  C/T3  resulting  from  the  different  crystallinities   are 
subscripted  to   distinquish  them. 

Examination  of   the  table  shows   that  the  value  deduced   from  this 
work  with  the  density  determined  crystal linity   is   in  excellent   agree- 
ment with  Isaacs   and  Garland's,    if    it    is   recognized  that   they  could 
not   separate  the  excess   heat   capacity.      Use  of  this   crystallinity 
rather  than  the  X-ray  value  contains   a  certain  amount  of  consistency 
since  the   present   results  were  based  on  the  density  measurements.      It 
is   noted   they  give  no  estimate  of  the  possible  error   in  the  X-ray 
value,    although    it    is  doubted  that    it   is  as  high  as   6%. 

Table  IV 

Comparison  of  present  results  with  those  of  Isaacs  and 
Garland.   (Values  expressed  in  mj/°lC-gm) 


C^r3(I&G) 

C/T3(90%) 

C/T3(85%) 

T  = 

=  2°K 

.0113 

.0097 

.0108 

3 

,0113 

.0100 

.0112 

4 

.0113 

.0102 

.0115 

5 

.0113 

.0103 

.0116 

(b)  Reese  and  Tucker  [  34  ] 

The  results  of  thermal  conductivity  measurements  on  three  sam- 
ples of  polyethylene  were  used  by  Reese  and  Tucker  to  indirectly 
determine  the  specific  heats  of  the  samples.   This  led  to  a  low  pre- 
cision averaged  set  of  values  over  the  temperature  range  1  to  4.5°K. 
The  comparison  with  the  present  work  was  made  in  the  same  manner  at 
it  was  for  Isaacs  and  Garland's  work,  and  equation  4.1  was  again  used, 


53 


However  in  order  to  compare  the  values,  a  conversion  was  necessary. 
All  of  the  results  of  Reese  and  Tucker  were  in  terms  of  volume  per- 
centages of  crystallinities  and  the  specific  heats  were  in  units  of 
mj/OK^-cmS.   The  values  used  below  are  not  as  reported,  but  have  been 
converted  using  their  values  of  the  densities. 

The  present  investigation  was  not  independent  of  the  work  of 
Reese  and  Tucker,  and,  in  fact,  was  an  outgrowth  of  it.  Therefore 
the  same  types  of  polyethylene  were  used  in  both  experiments.   How- 
ever slightly  different  densities  (crystallinities)  were  found  for 
the  samples  of  the  two  investigations.   This  is  explained  by  the  dif- 
ferent molding  processes,  and  cooling  times,  which  is  known  to  effect 
the  degree  of  crystallinity  [31  J. 

The  comparison  is  shown  in  Table  V.   The  sample  names  are  from 
the  work  of  Reese  and  Tucker  and  should  not  be  confuted  with  the  same 
names  used  in  the  present  work.   The  calculated  crystallinities  are 
shown  in  the  table. 

These  results  are  in  agreement  within  the  reported  error  given 
by  Reese  and  Tucker  (10%)  except  for  Marlex  I.   However  in  view  of  the 
indirect  methods  of  their  measurements,  it  is  not  considered  exces- 
sive.  In  addition  they  report  a  "dish  shape"  characteristic  if  the 
values  of  C/T  is  plotted  against  temperature.   This  was  proposed  to 
be  due  to  a  lack  of  precision  in  the  temperature  and  thermal  conductance 
determinations  at  the  two  temperature  extremes  of  their  range.   If  it 
is  assumed  that  this  was  the  case  at  only  the  low  end,  then  the  points 
above  2°K  will  display  one  half  of  the  "dish"  and  will  appear  as  the 
curves  in  figure  3  between  2  and  4.5°K.  This  is  believed  to  have 
been  the  case  since  at  the  lowest  temperatures,  the  parameter  from 
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which  the  heat   capacities   was   deduced   (the  thermal   equilibrium  times) 
would  be  greatly   increased  due  to   the  thermal   resistance   at  the  sam- 
ple-low temperature  sink   interface.      In  light   of   the  above,   the 
results   are  considered  to   be   in  good  agreement. 


Table  V 


Comparison  of  present  results  with  those  of  Reese  and 
Tucker  (Values  expressed  in  mj/°K  -gm) 


T  (°K) 


C/T3 


2 
3 
4 
5 

2 
3 

4 
5 

2 
3 

4 

5 


(Present) 


CA°(R&T) 
Marl  ex  I  (X  =  79%) 

.0147 
.0147 
.0147 
.0147 
Marl  ex  II  (X  =  62%) 
.0166 
.0166 
.0166 
.0166 
Low  Density  (X  =  47%) 
.0210 
.0210 
.0210 
.0210 
(c)  Sochava  [44] 

As  mentioned  in  the  introduction,  the  specific  heat  measure- 
ments reported  by  Sochava  for  the  temperature  range  17  to  60°K  are  on 
a  sample  of  unknown  density.   Since  it  has  been  shown  here  that  the 


.0118 
.0124 
.0128 
.0129 

.0149 
.0159 
.0167 
.0168 

.0175 
.0190 
.0200 
.0202 
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specific  heat  below  30°K  is  strongly  density  dependent,  comparison 
is  difficult.   However,  if  Sochava's  data  are  placed  on  a  specific 
heat-crystallinity  plot  from  the  present  work  and  are  adjusted  to 
fall  on  the  straight  lines,  a  crystallinity  of  53%  is  determined  for 
Sochava's  sample.   This  is  a  reasonable  value  for  commercial  poly- 
ethylene, and  it  is  felt  that,  as  far  as  the  comparison  can  be  made, 
there  is  agreement  between  the  two  sets  of  data, 
(d)  Dainton,  Evans,  Hoare,  and  Melia  [l3] 

Dainton,  et  al ,  reported  measurements  on  the  specific  heat 
of  two  samples  of  polyethylene  for  the  temperature  range  20  to  300°K. 
Their  two  samples  were  of  very  similar  crystallinity  (78  and  79%), 
and  due  to  the  manner  in  which  the  data  were  presented,  comparison 
could  be  made  only  at  20  and  30OR.   Again  the  crystallinity  was  de- 
duced from  the  reported  density  (it  was  given  for  only  one  sample), 
but  in  this  case  it  was  identical  with  that  given  as  determined  by 
X-ray  methods.   To  make  the  comparison  predictions  from  the  present 
results  had  to  be  adjusted  to  78  or  79%  crystall inities.   This  was 
accomplished  by  plotting  the  specific  heat  versus  crystallinity  for 
the  three  measured  samples  at  20  and  30°K,  connecting  the  points  with 
a  straight:  line,  and  finding  the  specific  heats  corresponding  to  the 
reported  crystallinities.   Table  VT  shows  the  results  of  this  predic- 
tion, along  with  the  results  of  Dainton' s  et  al .   The  subscripts  78% 
and  79%  refer  to  the  two  samples. 


56 


Table  VI 

Comparison  of  present  results  with  those  of  Dainton, 
et  al .   (Values  given  in  mj/°K-gm) 

Dainton,  et  al  Present  work 

c78%     c79%  c7  8%     c79% 

T  =  20°K      66.8     62.1  62.2     61.6 

30       147.0     146.2  146.7     146.4 

Considering  that  the  20°K  results  reported  by  Dainton,  et  al,  are 
extrapolated  from  22°K  (the  lowest  temperature  of  actual  measurements), 
the  agreement  is  felt  to  be  excellent. 

The  Tarasov  Model  [  47 ] 
The  present  results  show  clearly  that  a  strong  specific  heat- 
crystallinity  dependence  is  evident  in  polyethylene  below  30°K. 
Hence  no  single  set  of  Tarasov  parameters  can  be  expected  to  predict 
the  heat  capacity  in  general,  and  therefore  a  choice  was  made  to  sep- 
arate the  contributions  df  the  two  phases,  amorphous  and  crystalline, 
and  obtain  two  sets  of  thetas.   But  they  are  not  completely  indepen- 
dent, as  will  be  discussed. 

The  density  dependence  does  not  extend  much  beyond  the  tempera- 
ture range  of  this  investigation.   This  was  determined  from  the  fact 
that  the  results  of  Sochava  and  Dainton,  et  al ,  are  essentially  the 
same  above  50°K,  whereas  they  had  differed  by  over  10%  at  30°K.   Addi- 
tionally data  from  Dainton,  et  al ,  on  a  sample  of  crystallinity  58% 
above  90°K  show  that  the  specific  heat  for  this  sample  is  the  same  as 
for  the  more  crystalline  ones.   (Above  110°K  this  is  no  longer  true, 
but  it  is  due  to  effects  which  are  of  no  concern  here).   This  com- 
pares well  with  the  indications  of  figure  5  which  is  a  plot  of  the 
relative  difference  between  the  amorphous  and  crystalline  results 
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Figure  5.  Relative  difference  between  the  0%   and 
100$  extrapolated  data  versus  temperature. 
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below  30°K.   From  the  figure  it  is  clear  that  the  difference  will  be 
quite  small  for  temperatures  above  about  60°K,  if  the  curve  is  extra- 
polated above  30° K. 

The  disappearance  of  the  density  dependence  at  the  higher  temper- 
atures has  a  rather  natural  explanation  also  on  the  basis  of  the 
Tarasov  model.   The  three  dimensional  parameters  associated  with  the 
two  phases  are  small,  which  implies  that  the  three  dimensional  modes 
will  be  nearly  all  filled  at  not  very  high  temperatures.  Thus  the 
major  contribution  to  the  specific  heat  should  be  due  to  the  one 
dimensional  vibrations,  which  are  assumed  the  same  for  both  phases. 

The  assumption  that  the  one  dimensional  parameters  are  the  same 
in  both  the  amorphous  and  crystalline  cases  is  based  on  the  fact  that 
these  parameters  are  associated  only  with  the  forces  along  the  chains 
which  should  be  independent  of  the  configuration  the  chain  has  assumed. 
It  is  also  within  the  temperature  range  in  which  the  density  dependence 
disappears  that  the  specific  heat  becomes  linear  in  temperature,  as  is 
predicted  by  the  Tarasov  model,  and  it  was  from  this  temperature  re- 
gion that  the  one  dimensional  thetas  were  determined.   The  linear 
portion  of  the  specific  heat-temperature  curve  yielded  only  an  aver- 
age value  of  the  thetas,  and  more  information  was  needed  to  separate 
it  into  the  longitudinal  and  transverse  parts. 

The  assumed  value  of  1655°K  for  the  longitudinal  one  dimensional 
theta  was  based  on  the  measured  maximum  frequency  associated  with 
the  C-C  stretch  of  the  polyethylene  chain  as  determined  by  Raman 
spectrum  [ 15  J .   Since  this  parameter  had  been  taken  to  be  related  to 
the  maximum  acoustic  frequency  along  the  chain,  the  use  of  what  appears 
to  be  an  optical  mode  may  not  be  valid,  but  this  seemed  physically 
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representative  of  what  this  mode  should  be.   In  any  event  it  should  be 
pointed  out  that  the  predicted  specific  heat  on  the  extended  Tarasov 
model  is  relatively  insensitive  to  variations  of  the  values  of  &i~r 
and  "d-i«-  as  long  as  the  relation 


4.2 


r  -  "5"  (  %r   +   -&.L-) 


is  satisfied.   For  example,  if  the  ratio  of  "o>t  to  ~&7t-    is  taken  to 
be  .6  (the  same  ratio  as  assumed  for  the  three  dimensional  case),  it 
changes  the  specific  heat  preduction  at  20°K  by  only  .5%. 

One  further  remark  must  be  made  on  the  one  dimensional  para- 
meters.  On  the  basis  of  experimental  data,  obtained  by  neutron 
diffraction,  the  upper  limit  of  the  bond  stretching  acoustic  mode 
observed  in  polyethylene  was  found  to  be  545  t   25  cm-1  [ 32 J  -   This 
corresponds  to  a  characteristic  temperature  of  785  i  35°K.   If  this 
value  is  assumed  for  the  "#}  ,  it  predicts  specific  heats  which  are 
much  too  high  from  about  12°K  on  up  in  temperature.   This  is  inter- 
peted  not  as  anything  wrong  with  the  value  of  the  maximum  frequency 
noted  above,  but  rather  an  indication  that  there  is  no  real  physical 
significance  of  the  characteristic  temperature  from  the  continuum 
model . 

The  three  dimensional  thetas  for  the  amorphous  and  crystalline 
parts  differ  from  one  another  in  the  manner  in  which  they  should  be 
expected.   They  are  related  to  the  forces  acting  between  the  poly- 
ethylene chains,  and,  since  the  forces  are  most  probably  of  the  type 
which  decrease  rapidly  with  distance,  the  less  dense  and  less  closely 
packed  chains  of  the  amorphous  state  should  have  smaller  thetas.   The 
three  dimensional  parameter  deduced  from  the  experimental  data  support 
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this  qualitative  description. 

The  calculated  values  of  the  specific  heat  based  on  the  original 
Tarasov  model  are  not  in  good  agreement  with  the  experimental  re- 
sults.  As  can  be  seen  in  Appendix  D,  there  is  a  5-6%  difference 
around  20°K,  although  in  other  temperature  ranges  the  fit  is  better. 
The  extended  Tarasov  model,  which  separated  the  longitudinal  and 
transverse  contributions,  was  used  since  it  gave  an  approximately 
2-4%  improvement  over  the  entire  range  when  compared  with  the  basic 
model.   However  in  the  final  analysis  the  Tarasov  model  should  not 
be  expected  to  explain  the  specific  heat  of  polyethylene  any  more 
than  the  Debye  theory  can  for  other  solids,  since  they  contain  inher- 
ent erroiieous  assumptions  concerning  the  lack  of  dispersion.   In  view 
of  this,  the  fact  that  the  Tarasov  model  fits  as  well  as  it  does 
makes  it  useful  as  a  first  approximation  for  the  specific  heats  of 
fibrous  solids  such  as  polyethylene  at  low  temperatures. 
The  Stockmayer-Hecht  Model  [ 19,46  ] 

In  order  to  investigate  the  applicability  of  the  model  of  Stock- 
mayer  and  Hecht  as  expanded  by  Genensky  and  Newell,  the  specific  heat 
of  the  100%  crystallinity  case  was  the  only  one  considered,  as  the 
model  was  suggested  for  a  fibrous  crystalline  solid.   At  the  lowest 
temperatures,  the  model  did  predict  the  cubic  temperature  dependence, 
but  the  predicted  temperature  region  for  which  the  specific  heat  would 
be  proportional  to  T  *   was  not  found.   However  it  is  felt  that  this 
was  an  over  simplification  on  the  part  of  Genensky  and  Newell  for, 
even  if  their  analysis  were  correct,  it  would  not  be  rigorously  re- 
flected by  this  relation.   The  reason  for  this  is  that  they  predict 
in  some  temperature  region  the  specific  heat  should  be 
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4.3 


c  -.    Ht      -+  &r' 


where  A  and  B  are  constants  (the  first  term  is  the  contribution  of 
the  interchain  interactions  and  the  second  is  the  result  of  vibrations 
along  the  chains).   This  prediction  does  not  give  any  appreciable 
temperature  region  in  which  the  specific  heat  is  proportional  to  a 
constant  power  of  the  temperature. 

To  test  the  model  for  the  temperatures  at  which  the  specific 
heat  was  investigated,  a  plot  of  C/T  '   versus  the  square  root  of  the 
temperature  should  have  a  region  where  the  slope  is  constant,  if  the 
model  is  correct.   Figure  6  shows  that  this  does  not  occur  in  the 
temperature  range  below  30°K,  and,  since  the  slope  is  negative  at 
the  high  temperature  end,  it  is  doubted  that  it  will  occur. 

However,  as  is  also  shown  in  figure  6,  plotting  C/T  versus  T  ' 
does  show  a  linear  region  (22  to  27°K)  which  imply s 

5-/2. 


4.4 


C  =  AT      +ST. 


This  could  be  interpeted  as  implying  that  their  analysis  is  incorrect 
in  the  sense  that  they  do  not  anticipate  that  the  linear  temperature 

contribution  from  the  vibrations  along  the  chain  would  occur  at  the 

5/2 
same  temperatures  that  the  T    relation  holds  for  the  interchain 

specific  heat  part.   However,  this  cannot  be  reconciled  with  the  model 

as  presented,  because  the  relative  size  of  the  two  contributions  is 

wrong. 

It  is  not  possible  in  practice  to  draw  from  experimental  heat 

capacity  data  reliable  inferences  regarding  the  detailed  form  of  the 

frequency  distribution,  as  the  heat  capacity  tends  to  smooth  out  the 

individual  features  of  the  frequency  distribution.   Hence  the  compar- 
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ison  with  the  model  of  Stockmayer  and  Hecht  cannot  specifically  deter- 
mine where  their  analysis  is  in  error.   On  the  basis  of  the  apparent 
T  /  relation  for  the  interchain  specific  heat  contribution  (if  the 
interpetation  of  figure  6  is  correct)  ,  it  seems  that  the  essence  of 
the  model  is  correct,  but  that  a  reexamination  of  the  one  dimensional 
part  is  in  order.   However  other  assumptions  of  theirs  may  be  wrong. 
Since  the  polyethylene  crystal  appears  to  be  lamellae,  (although  not 
directly  observed  in  the  bulk  polyethylene)  and  the  crystalline  form 
is  spherulitic,  the  structure  may  be  such  that  the  idealized  tetra- 
gonal crystal,  assumed  by  Stockmayer  and  Hecht,  will  not  yield  the 
correct  predictions.   Moreover  since  they  had  to  make  estimates  of 
the  relative  sizes  of  the  force  constants,  of  which  little  is  known, 
erroneous  assumptions  here  could  lead  to  different  forms  for  the  fre- 
quency spectrum,  and  thereby  errors  in  the  deduced  heat  capacity.   In 
any  event  in  its  present  form  the  specific  heat  predicted  from  this 
model  is  not  completely  applicable  to  polyethylene  below  30°K. 

The  Excess  Heat  Capacity 
In  the  absence  of  any  acoustical  measurements  made  at  low  temp- 
eratures for  the  heat  capacity  of  polyethylene,  the  humping  of  the 
C/T  versus  temperature  curve,  figure  4,  was  interpeted  as  an  excess 
heat  capacity.   The  appearance  of  this  excess  was  not  expected,  but 
it  is  not  unreasonable  that  it  does  appear  when  it  is  recognized  that 
it  is  present  in  other  amorphous  solids  [4,12,17J.   The  origin  of  the 
excess  heat  capacity  i*  unknown,  but  attributing  it  to  the  existence 
of  non-acoustical  low  frequency  modes  seems  reasonable.   In  particu- 
lar, the  Rosenstock  theory  is  appealing  since  it  is  known  that  free 
volume  (holes)  does  exist  in  amorphous  polyethylene  [  10, II*].   However 
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no  reasonable  conclusions  can  be  reached  on  this  or  any  other  model 
which  might  explain  the  excess  on  the  basis  of  the  present  results. 

The  addition  of  the  Einstein  term  to  the  amorphous  heat  capacity 
empirically  fits  the  data  quite  well  at  the  low  temperatures,  and  no 
other  correction  is  deemed  necessary  to  account  for  the  excess  in  an 
analytic  expression  for  the  specific  heat.   In  fact,  it  is  felt  that 
when  low  temperature  acoustic  measurements  of  the  specific  heat  of 
polyethylene  become  available,  the  value  for  the  extrapolated  T  =  0°K 
point  will  agree  with  the  following  relation  which  is  based  on  the 
present  results; 

From  the  size  of  the  constant  which  multiplies  the  Einstein 
function,  it  is  readily  deduced  that  of  all  of  the  possible  vibration 
modes  of  the  solid,  only  . 17%  are  taken  up  in  the  excess  term.   This 
can  be  compared  with  the  number  of  these  modes  in  vitreous  germanium 
(.4%)  and  vitreous  silica  (2%),  but  no  conclusions  can  be  drawn  [4,17], 

Summary 

As  has  been  shown,    both   the   continuum,    Tarasov,  model   and  the 
lattice  dynamical,    Stockmayer-Hecht ,   model   fail   to   completely  predict 
the  specific  heat  of   polyethylene.      This   is   not   surprising   as  they 
are  based  on  models   which  cannot  be  exact   for  differing  reasons.      The 
continuum  model    is  at   best   a  first   approximation   since   it   ignores   the 
dispersion  of  the   sound  velocity  which   surely  takes   place   in  any  real 
solid.      Hence  the   spectrum  of  vibrational   states    is  wrong.      But,    con- 
sidered  as   a  simple   two   parameter   approximation,    the  predictions  are 
very  good.      The  Stockmayer-Hecht  model   which  attempts   to   account   for 
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the  dispersion  fails  for  different  reasons.   Since  this  model  is 
based  on  an  assumption  of  the  interparticle  forces  in  various  direc- 
tions, there  are  many  parameters  which  determine  the  predictions.   If 
some  are  chosen  incorrectly,  the  whole  model  fails.   On  the  basis  of 
the  present  results  no  attempt  was  made  to  evaluate  the  parameters, 
since  one  of  the  basic  predictions  of  the  model  did  not  materialize. 

The  Einstein  term  adjustment  for  the  excess  heat  capacity  found 
in  the  amorphous  polyethylene  was  adequate  to  explain  the  departure 
from  the  cubic  temperature  dependence  of  the  specific  heat  at  the 
lowest  temperatures.   But  until  such  time  that  empirical  evidence  is 
brought  forward  to  verify  some  model  to  explain  it  completely,  it  must 
remain  as  an  a  posteriori  addition  to  explain  the  experimental  re- 
sults. 

The  hypothesis  that  the  specific  heat  of  polyethylene  is  linear 
in  crystallinity  was  proven,  it  is  felt,  in  the  temperature  range 
covered  by  this  investigation,  and  extrapolation  to  the  completely 
amorphous  and  completely  crystalline  cases  is  considered  valid.   It 
is  also  believed  that  the  crystallinity  dependence  has  disappeared  by 
50-60°K  (considering  only  the  temperatures  below  110°K). 

The  present  low  temperature  specific  heat  measurements  are  more 
reliable  than  those  given  by  the  prior  investigations  (Isaacs  and 
Garland,  and  Reese  and  Tucker),  since  the  existence  of  the  excess 
heat  capacity  is  recognized.   By  extrapolating  these  results  to  0°K 
and  by  using  the  method  of  interpolating  between  crystallinities 
noted  before,  this  work  can  be  mated  to  that  of  Dainton,  et  al ,  to 
give  a  continuous  specific  heat-temperature  coverage  up  to  HO°K, 
where  many  other  data  become  available.   The  resulting  information 
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should  prove  valuable  to  workers    investigating  models  of  crystalline 
and   non-crystalline  fibrous   solids    into  which  nearly  all   polymers  can 
be  classed. 
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APPENDIX  A-  Calibration  of  the  Thermometers 
Introduction 

Probably  the  most  exacting  aspect  of  this  experiment  was  the 
determination  of  the  absolute  temperature  scale.   The  difficulty 
stems  from  the  fact  that  between  5  and  14°K  there  is  no  convenient 
temperature  standard.   Below  5°K  the  scale  had  been  set  by  inter- 
national convention  in  terms  of  the  vapor  pressure  of  liquid  helium, 
and  between  14  and  25°K,  the  vapor  pressure  of  liquid  hydrogen  is  a 
good  standard  [8,48,51].   Additionally  the  internationally  accepted 
platinum  resistance  thermometers  can  be  used  above  10°K  [20J.   Be- 
tween 5  and  14°K  the  temperature  scale  is  normally  established  by 
means  of  a  gas  bulb  thermometer. 

In  this  investigation  the  resistance  thermometers  were  calibra- 
ted below  4.2°K  by  means  of  a  vapor  bulb  thermometer,  and  from  4.2 
to  30°K  by  means  of  a  gas  bulb  thermometer.   Liquid  hydrogen  vapor 
pressure  thermometery  was  not  attempted  due  to  its  inherent  dangers: 
the  facilities  were  not  available  to  obviate  them.   Initially  it  was 
planned  that  the  platinum  resistance  scale  was  to  be  compared  with 
the  gas  bulb  results,  but  was  not  since  another  comparison  became 
available  (the  specific  heat  of  copper,  Appendix  B),  and  it  soon  ap- 
peared that  the  higher  accuracy  in  temperature  range  above  10°K  was 
not  necessary  for  the  purpose  of  this  investigation. 

The  theory,  construction  and  experimental  procedures  concerned 
with  the  calibration  of  the  resistance  thermometers  is  the  subject  of 
the  remainder  of  this  appendix,  together  with  an  estimate  of  the 
error  in  the  temperature  determination. 
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Theory  of  Constant  Volume  Gas  Bulb  Thermometry 
The  equation  of  state  for  an  ideal  gas  is, 

a.i  «£«■&& 

Ti  T+ 

where  a  and  /  refer  to  the  initial  and  final  states  of  the  system. 
If  the  volume  is  held  constant  and  the  initial  conditions  are  known, 
then  there  exists  a  linear  relation  between  the  temperature  of  the 
gas  and  its  pressure.   A  constant  volume  gas  bulb  thermometer  could 
consist  of  a  large  volume,  V  ,  at  the  temperature  to  be  measured,  '  , 
and  a  pressure  sensing  device  at  temperature,  'y   ,  with  volume  ^con- 
nected by  a  tube  of  small  volume  ^*  .   Then  equation  A.I  must  be  modi- 
fied to  account  for  the  fact  that  not  all  the  gas  is  at  the  low 
temperature  in  the  final  condition;  some  is  still  in  the  gauge  and 
connecting  tube.   To  account  for  this  equation  A.I  is  replaced  by, 


<i*T 

where  i  y  is  the  pressure  of  the  system  when  it  is  at  Ty  ,  and  i  is 

the  pressure  when  the  bulb  is  at  temperature  '  .   The  integral  ex- 
presses the  fact  that  the  temperature  along  the  connecting  tube 
varies  from  '  to  ly  .      However,  it  is  convenient  to  assign  a  temper- 
ature Tc.  such  that, 

)       J.  f  oitr'Ct) 

A'3        ~TiS     *r'    I  T~    * 

Solving  equation  A. 2   for "T ,    it   is   found  that, 

a.*    t  *  tv  p  -—-, — y. 
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Consider  now  the  case  where  If^O,    that   is  assuming  that  the 
connecting  tube's  volume    is   quite  small    in  comparison  with  that  of 
the  other  components.      Then  equation  A.  4  becomes, 

en 

A. 5  T  = ; 7Z     ' 

Then  if  -^"Vs  /  ,    and    if  the  bulb   is   at    low  temperature,   most  of  the 
gas  will  be  in  bulb,    and   the  pressure  will  be  essentially   linear    in 
temperature. 

However,    in  an  actual   experiment   the  gas  will   not  be   ideal,    as 
implied   in  the  previous  arguments.      The  equation  of    state   in  this 
case  can  be  written, 

where  the  8< *   are  the  virial   coefficients,   and  LT  is   the  specific 
volume.      If  helium  were  used    in  the  temperature  range  of  the   present 
experiment   8>  is  of  2  or   3  orders  of  magnitude  smaller  than   8/    ,    and 
the  term  containing    Bz  can  be  neglected  [  22  ].      Since  8,  is   also   small, 
t^can  be  set   equal  to  ~~p~  in  the  term  which  corrects   for  the  second 
virial   coefficient,    so  that, 

A.7         f>V  =  *l?(-Ti-Bf) 


putting  -=sr  "  &  • 

Using  A.7  equation  A. 2  becomes, 


p  (      v  i       "~      -*      °-        1  ^  P      v       J-  £      + 


i  | 

where  the  subscripts  on  the  S>  f  refer  to  the  temperature  at  which 
they  are  evaluated. 

Finally,  consider  the  form  that  equation  A. 8  would  take  if  the 
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initial  conditions  were  varied  somewhat,  that  is,  by  filling  the 
system  with  the  gas  when  the  temperature  of  the  bulb  is  at  a  third 
temperature  ~TN   and  therefore  has  volume  r*    due  to  thermal  contrac- 
tion.  Under  these  circumstances,  equation  A. 8  will  be, 


p(   v  -  -+  — £—    -i  ^lL 

\T-tBP  ^  t,+b,p  n  rc'+8c-p 

where  the   fv  refers   to   the   pressure   in  the   system  when  the  bulb   is   at 
temperature  77,  •      (Notice  that   the   Tc  s  on  the  two   sides  of  the  equa- 
tion are  not  the  same  since  the   lower  temperature   is   different.) 
Construction  Gas   Bulb  Thermometer 

As    indicated   in  the  comment   following  equation  It. 5,   temperature- 
pressure   linearity  can  be  obtained    if  the  gauge  volume  to   bulb  volume 
ratio    is  made  small.      As   a  design  criteria   -rr-    was   taken  as   1/20,    and 
the  connecting  tube's  volume   IT   would  be  held   small  compared  to  the 
volume  of   the  gauge. 

To   obtain  the  desired  temperature  range,    the  bulb  was    initially 
cooled  to   4.2°K  by  a   liquid  helium  bath,    and  then  heated  electrically 
above  this   point.      The  bulb  was  made  of  copper,    to   insure   uniformity 
of  temperature  throughout    it,    and  the  resistors  were  placed  directly 
on  the  outer   surface  of   the  bulb. 

The  pressure  measuring  device  was  a  Wallace  and  Tierman  gauge, 
type  FA-145    (2   turn  dial,    0-800  mm.    Hg   range).      A    .44  mm.    inside 
diameter   stainless   steel  capillary  tube  served   as   the  connecting 
tube,    and   the  gas   used  was  helium. 

The  gauge  was   calibrated   against   a  mercury  manometer    (assumed 
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correct)  both  before  and  after  calibration  of  the  resistors.   There 
appeared  a  linear  deviation  with  pressure  which  was  regarded  as  a 
systematic  error  in  the  gauge  and  a  correction  for  it  was  made  in  the 
subsequent  calculations. 

The  volume  of  the  gauge,  ^,  was  determined  by  comparison  with  a 
known  volume.   The  value  was 

IT   -  10.4(1  +  1.4  X  10-4P)  1  .3  cm3 
where  P  is  the  pressure  in  mm.  of  Hg.   (This  correction  term  accounts 
for  the  fact  that  the  gauge  had  a  bellows  type  pressure  sensing  ele- 
ment) . 

Since  it  was  desired  that  V      the  volume  of  the  bulb  should  be 
at  least  230  cm3.   A  1  1/2"  copper  pipe,  10"  long,  sealed  at  the  ends 
produced  a  volume  calculated  to  be  273.4  t    .7  cm3  at  room  temperature. 

Because  of  the  size  of  the  dewar,  it  was  necessary  to  have  the 
capillary  48"  long.  This  resulted  in  a  calculated  volume  for  the 
connecting  tube  of  only  .19  cm3,  well  within  the  design  criteria. 

As  liquid  helium  temperature  (approximately  4.2°K)  would  be  the 
lowest  temperature  used  in  the  calibration,  some  means  must  be  de- 
vised to  thermally  isolate  the  bulb  from  the  bath,  once  it  was  in 
equilibrium  with  the  bath,  in  order  to  raise  it  above  this  point. 
For  this  reason  the  bulb  was  suspended  in  a  vacuum  container  which 
was  to  be  immersed  in  the  bath.   The  bulb  support  was  a  3/8"  O.D. 
(.01"  wall  thickness)  stainless  steel  tube.   This  tube  also  extended 
up  and  out  of  the  cryostat  and  served  as  the  pumping  line.   An  elect- 
rical heater  of  about  1500  ohms,  consisting  of  12  feet  of  .002"  man- 
ganin  wire  (f  =  78.3  ohms/ft.)  was  wound  around  the  bulb  and  glued 
there  with  G.  E.  7031  varnish.   The  heater  was  necessary  to  raise  the 
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temperature  and  to  compensate  for  the  resultant  heat  leak  through  the 
support  (and  to  a  much  lesser  extent,  the  heat  leak  due  to  radiation 
and  conduction  by  any  residual  gas). 

The  resistors  to  be  calibrated  were  inserted  into  small  recep- 
tacles in  the  bulb.   In  order  to  avoid  the  relatively  large  electrical 
lead  resistance,  a  three  wire  connection  was  made  to  each  resistor. 
This,  and  a  discussion  of  the  low  power  A.  C.  bridge  used,  is  given  in 
the  work  by  Reese  and  Tucker,  and  will  not  be  presented  here  [34]. 

A  schematic  diagram  of  the  completed  system  in  shown  in  figure  7. 
Gas  Bulb  Calibration  Procedures 

The  gas  bulb  was  filled  with  gaseous  helium  to  about  790  mm.  of 
Hg  when  the  system  was  at  liquid  air  temperature.   The  vacuum  spaces 
were  then  pump#d  to  approximately  10""'  Torr  and  the  liquid  helium 
space  was  filled. 

The  data  accrual  was  quite  simple.   The  procedure  was  to  heat 
for  a  period  of  about  30  seconds,  then  turn  off  the  heating  power,  and 
then  observe  the  resistances  come  to  equilibrium,  which  indicated 
temperature  equilibrium.   The  pressures  and  resistance  values  were 
then  recorded.   The  tendency  for  the  temperature  to  drift  downward, 
due  to  the  heat  leak,  was  compensated  for  by  the  heater. 

The  observations  were  taken  over  a  pressure  range  of  about  40- 
300  mm.  of  Hg,  corresponding  to  temperatures  of  from  4  to  30  degrees. 
They  were  taken  twice  over  the  entire  range,  at  about  2-4  mm.  steps, 
to  check  reproductibility ,  and  a  third  time  over  the  lowest  portion 
of  the  range  where  the  resistors  were  the  most  sensitive.   Seventy- 
three  points  were  taken  altogether. 

After  completion  of  this  procedure,  the  bulb  was  allowed  to  come 
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To  manometer 


/       To  vacuum  sources 


a.  Vapor  bulb 

b.  Resistors  to  be  calibrated 

Figure  7.   Schematic  of  gas  bulb  thermometer. 
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to  complete  equilibrium  with  the  liquid  helium  bath.   Using  the  vapor 
pressure,  the  bath  temperature  was  determined.   This  would  later 
serve  as  a  check  on  the  accuracy  of  the  calibrating  procedures.   The 
helium  was  then  evaporated  from  the  dewar  and  the  space  filled  with 
liquid  nitrogen  (99.9%  pure).   The  pressure  in  the  gas  bulb  together 
with  the  known  temperature  of  the  nitrogen  bath  gave  values  of  f* 
and  T*  which  could  then  be  used  in  equation  A. 9,  and  the  pressures 
could  be  converted  into  temperatures. 

Equation  A. 9  is  here  rewritten  with  a  slight  amplification  of 
the  terms. 


f»  Vm  ,      ft.  tr  ,  P-  v ' 


A. 10 


■f 


T-tfSP         Tyii3yP       -rj+£? 

where   Y»    is   the  pressure   in  the   system  when  the  bulb   is  at   liquid 
nitrogen  temperature,   T^   is   the  liquid   nitrogen  temperature,    V     is 
the  volume  of   the  bulb  at   temperatures   less  than   30°K,    and    V#   is  the 
volume  of  the  bulb  at  T*  .      The  subscript  r   refers   to  ambient   condi- 
tions in  the   laboratory,    and  the  other   symbols  have  the  same  meaning 
as   before.      Using  published  values   for  the  thermal   contraction  of  cop- 
per at  the  various  temperatures,    it   is   found  that   Vr 270.8 t  .7   cm3  and 
Vv=  273.0  t  .7   cm3   [26].      For  an  estimate  of  the   "^ 's  ,    a  value  of   145°K 
was   assumed   foT    'c    ,   and   108°K  for  "FT  .      These  were  obtained  by  assuming 
a   linear  temperature  distribution  along   the  capillary.      This    is   suf- 
ficiently accurate  as  this   part  of  the  system   is  quite  small.      Values 
of  the  second  virial   coefficients  were   taken  from  Keeson;(_22j 
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A. 11  0    =    (.2453  -    6.2009/T)    °K/mm  of  Hg 

fy  =    .000163 

8y  =    .000162 

Qc  =    .000179 

6C  =.000186    . 
Using  the  above   in  equation  A. 10   leads   to   an  equation  quadratic   in 
Solving  that    for    T yields, 


where 


A.i2      t=  --£  +  i- /*"-*; 

-*-(£- -Op 

Equation  A. 12  was  programmed   for   the  CDC   1604  digital   computer. 

The  procedure  discussed   above  gave  only  a  pressure-resistance 
relationship  until   equation  A. 12  was   solved.      To  obtain  the  desired 
temperature-resistance  dependence,    all   the  data  was   best   fitted,    in 
the   least   squares   sense,    by  the  computer  to  the   following  type  of 
equation  for  each  resistor; 

3,      ,  ./*-;-< 

A. 13  T  =    £  A:  [JLy  R) 

where  the  A*  $  and  the    best  fitted  /*  (integers  between  1  and  7) 
were  determined.   The  values  are  presented  at  the  end  of  this  appen- 
dix, and  the  RMS  deviation  from  this  fitting  is  calculated  to  be  .23%. 
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Calibration  of  the  Thermometers  below  4.2°K 

For  temperature  determination  from  1  to  5.2°K,  the  International 
Commission  on  Weights  and  Measures  adapted  the  "1958  He^  Vapor  Pres- 
sure Scale  of  Temperatures"  as  the  international  standard  [8,48]. 
Tables  based  on  this  scale  give  a  conversion  from  vapor  pressures  to 
temperatures,  hence  the  temperature  can  be  ascertained  from  a  reading 
of  the  vapor  pressure  alone. 

For  this  lower  temperature  determination,  it  was  decided  to  use 
a  vapor  bulb  thermometer,  rather  than  merely  measuring  the  pressure 
above  the  bath  since  the  latter  procedure  involves  making  corrections 
and  assumptions  of  questionable  validity. 

Vapor  bulb  thermometery  is  performed  by  placing  a  container  into 
the  helium  bath  at  the  liquid  level  at  which  the  temperature  is  de- 
sired to  be  known,  and  filling  the  container  with  gaseous  helium 
until  the  gas  condenses  which  indicates  equilibrium  between  the 
liquid  and  gas  phases  at  this  pressure  [42],   The  pressure  is  then 
observed,  and  tables  give  the  temperature  directly. 

For  the  actual  construction  of  the  bulb  a  small  copper  cylinder 
with  a  volume  of  approximately  2  cmJ  was  connected  by  a  1/8'  stainless 
tube  to  a  mercury  manometer  and  filling  manifold.   The  bulb  was 
mounted  directly  on  top  of  the  vacuum  container,  and  connecting  tube 
was  brought  outside  of  the  cryostat  enclosed  in  an  evacuated  tube  (to 
prevent  gas  oscillations  in  the  system  [42].   See  figure  7. 

For  the  calibration,  exchange  gas  was  introduced  into  the  vacuum 
space  where  the  resistors  were  located,  and  the  pressure  above  the 
bath  was  lowered  in  steps.   The  manometer  readings  were  taken  and 
the  corresponding  resistances  were  recorded.   In  all  32  points  were 
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taken  in  the  temperature  range  of  from  2.4  to  4.2°K. 

After  conversion  of  the  pressures  into  temperatures,  the  resis- 
tance values  were  related  to  the  temperature  by  equation  A. 13. 
Recalibration  of  the  Thermometers 

Inasmuch  as  the  comparison  experiment  (Appendix  B)  indicated  that 
a  systematic  error  might  be  present  in  the  temperature  scale,  the 
scale  was  reestablished.   Also,  in  the  initial  calibrations,  the  vapor 
pressure  and  the  gas  bulb  temperature  scales  (below  and  above  4.2°K) 
were  independently  determined  with  the  result  that  there  existed  a 
slight  discrepancy  at  their  common  temperatures.   (Approximately  .2% 
difference).   Since  it  was  desired  to  have  a  continuous  scale  over  the 
entire  range,  the  calibration  was  repeated  insuring  that  the  tempera- 
ture overlap  was  consistant  with  both  types  of  temperature  determina- 
tion.  This  was  accomplished  by  using  the  liquid  helium  temperatures, 
as  determined  by  vapor  pressures,  as  the  fixed  points  for  the  gas  bulb 
calibration  rather  than  the  liquid  nitrogen  point.   As  a  check  as  to 
the  precision  of  this  method,  the  predicted  liquid  nitrogen  tempera- 
ture was  only  . 3°K  different  from  the  assumed  value  of  77.35°K,  a 
difference  of  only  .23%. 

The  equipment  and  procedures  were  the  same  as  outlined  earlier. 
A  total  of  99  points  were  taken  over  the  entire  range,  varying  the 
intervals  so  that  a  uniform  coverage  was  obtained.   The  gauge  volume 
was  again  experimentally  determined,  and  was  within  3%  of  the  original 
measurement.   The  gauge  was  likewise  recalibrated  against  a  mercury 
manometer  twice. 

The  resistance  values  were  related  to  the  temperature  by  use  of 
equation  A.  13. 
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Accuracy  of  the  Calibration 

In  the  temperature  range  below  4.2°K  where  the  vapor  bulb  ther- 
mometer was  employed,  the  only  significant  source  of  error  in  the 
temperature  determination  was  the  reading  of  the  manometer.   It  could 
be  read  to  within  .001  cm  of  Hg  with  reproducibility  of  .003  cm.  of 
Hg.   This  results  in  a  maximum  uncertainty  in  the  temperature  for  this 
range  of  3  millidegrees . 

For  the  gas  bulb  calibration,  which  covered  the  temperature  range 
of  from  4.2  to  30°K,  the  following  effects  were  considered  to  ascer- 
tain an  estimate  of  the  error; 

(a)  The  fixed  temperature  point. 

The  calibration  over  the  entire  range  is  very  sensitive  to 
any  errors  in  the  determination  of  the  fixed  temperature  point  which 
essentially  sets  a  value  for  the  right  hand  side  of  equation  A. 12.   To 
determine  this  parameter,  six  independent  observations  of  the  tempera- 
ture were  made,  by  vapor  bulb  methods,  at  the  liquid  helium  bath  tem- 
perature, and  the  pressure  associated  with  each  recorded.   The  fixed 
point  was  taken  to  be  the  average  of  these  data.   The  resulting  error 
in  the  temperature  determination  over  the  entire  range  due  to  this 
was  estimated  to  be  at  most  .00  3  degrees. 

(b)  Errors  in  the  measurement  of  volumes. 

Since  the  bulb's  volume  was  much  greater  than  the  other  vol- 
umes in  the  system,  the  temperature  is  essentially  linear  in  the 
pressure.   This  tends  to  make  any  error  in  the  measurements  of  the 
volumes  self  compensating,  and  their  net  contribution  to  the  error  in 
the  calibration  is  negligible. 
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(c)  Thermalmolecular  correction. 

Following  Roberts  and  Sydoriak,  the  thermalmolecular  cor- 
rection is  .004°  at  the  low  end  of  the  temperature  scale  (4.2°K), 
and  negligible  above  7°K  [37].   This  systematic  error  could  have  been 
corrected  for,  but  was  not  in  view  of  the  existing  larger  errors  which 
will  be  discussed  below. 

(d)  Pressure  measurement. 

The  pressure  measuring  gauge  readings  were  compared  with  those 
of  a  mercury  manometer  both  before  and  after  the  resistor  calibration. 
This  comparison  was  made  stepwise  going  up  in  pressures,  since  this 
was  the  way  data  were  taken,  to  include  any  hysteresis  of  the  gauge's 
bellows  as  part  of  the  correction.   The  gauge  was  determined  to  have 
a  linear  error  as  a  function  of  pressure,  that  is,  the  deviation  was 
consistently  high  at  the  higher  pressures,  and  low  at  the  lower  ones. 
The  data  were  corrected  for  the  presence  of  this  systematic  error. 
The  statistical  uncertainty  of  this  correction  is  entirely  compatible 
with  the  precision  to  which  the  gauge  could  be  read,  about  .1  mm.   Due 
to  the  uncertainty  in  the  pressure  determination,  a  possible  error  of 
.01°  exists  throughout  the  temperature  range. 

Other  sources  of  error,  such  as  those  considered  by  Beattie,  and 
those  due  to  the  uncertainty  of  the  values  of  the  virial  coefficients, 
accuracy  of  the  mercury  manometer,  etc.,  were  considered  to  be  insig- 
nificant compared  with  the  ones  discussed  above  L'5]. 

Once  the  temperature  scale  was  established,  there  exists  the 
necessity  of  relating  the  resistance  values  of  the  resistors  with  the 
temperatures.   This  introduces  another  possible  reading  error  which 
amounts  to  .  17o  over  the  entire  range.   However  it  is  felt  that  due  to 
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the  insensitivity  of  the  resistors  at  the  higher  temperatures,  a  more 
realistic  value  of  .3%  should  be  assigned  there.   After  the  resistance- 
temperature  relation  had  been  established  by  equation  A. 13,  it  was  pos- 
sible to  assign  the  maximum  temperature  differences  these  errors  would 
introduce,  .0015°  at  2.5  and  .03°  at  the  highest  temperatures.   It  is 
also  noted  that  to  use  the  resistors  as  thermometer,  another  reading 
error  of  this  size  could  exist. 

Summarizing  all  of  the  above,  and  recognizing  that  the  fitted 
relation  between  the  temperature  and  the  sample  thermometer  resis- 
tance values  contains  an  uncertainty,  it  is  felt  that  the  temperature 
determination  by  use  of  the  resistors  was  accurate  to  within  .5%  of 
the  absolute  temperature. 

As  an  estimate  of  the  reproducibility  and  precision  of  the  cali- 
bration, the  RMS  deviation  of  the  actual  points  from  the  fitted 
resistance-temperature  was  used.   This  amounted  to  .2  3%,  and  appeared 
to  be  randomized.   Hence  the  statistical  scatter  of  the  calibration 
points  from  the  fitted  relation  was  well  within  the  estimated  errors 
of  the  procedure. 
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Parameters   for  the  temperature-resistance  fitted  relation 


7=  I  A; (Ay**)' 


Thermometer 

Shield 

Shield 

Sample 

(4.2-30°K) 

(2.4-4.2QK) 

(2.4-30°K) 

n 

6 

6 

4 

Al 

.81986X10-1 

-.11574X10"1 

-.24787X10° 

A2 

.65400X10-2 

.66904X10-1 

.32888X101 

A3 

-.50749X101 

.25718X10-1 

-.10029X102 

A4 

.17961X102 

.10939X10° 

-.36394X102 

A5 

.14366X102 

-.27080X101 

.25491X103 

A6 

-.15249X103 

-.69149X101 

-.45361X103 

Ay 

.19880X103 

.38145X102 

.36772X103 

%  RMS 

Deviation     .20 

.00 

.23 
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APPENDIX  B  -  Comparison  Experiment 
There  has  been  considerable  discussion  about  the  need  for,  and 
choice  of,  a  suitable  calorimetric  standard  in  the  temperature  range 
below  20°K  [53].   Because  of  the  large  number  of  specific  heat  deter- 
minations for  copper,  its  available  purity,  and  other  reasons,  there 
is  much  support  for  its  choice  as  such  a  standard.   At  the  same  time 
that  this  experiment  was  being  performed,  Martin  reported  the  specific 
heat  of  a  cast  99.999%  pure  copper  sample  in  the  temperature  range  of 
this  experiment  [30  J.   He  asserts  that  his  results  are  accurate  to 
better  than  1%,  which  implys  a  temperature  determination  to  within 
. 35%  accuracy. 

To  ascertain  an  estimate  on  the  accuracy  of  the  temperature 
scale  used  in  the  present  experiment,  it  was  decided  to  duplicate 
Martin's  measurements.   A  single  crystal  copper  sample  of  99.9995% 
purity  was  obtained  (Semi-Elements,  Inc.),  and  was  etched  with  nitric 
acid  to  a  size  which  was  compatible  with  the  calorimeter.   The  re- 
sulting sample  was  a  2V  by  9/16"  rod  with  a  mass  of  83.181  grams. 
Heat  capacity  measurements  were  performed  on  this  sample,  in  the  same 
manner  as  for  the  polyethylene  samples,  over  the  temperature  range 
4.5  to  30°K. 

In  performing  these  measurements,  the  only  experimental  procedu- 
ral differences  from  the  polyethylene  runs  were  that  the  amount  of 
heat  introduced  for  the  temperature  steps  was  considerably  less,  and 
the  equilibrium  times  were  much  shorter,  due  to  the  differing  heat 
capacity  and  thermal  conductance.  Three  complete  sets  of  measure- 
ments were  obtained  over  the  aforementioned  temperature  range,  re- 
sulting in  a  total  of  119  data  points. 
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After  smoothing  the  data  by  least  squares  fitting  it  to  a  seven 
parameter  equation  (equation  3.2)  and  correcting  for  the  heat  capacity 
of  the  addenda,  the  results  were  compared  with  those  of  Martin.   They 
matched  fairly  well  at  the  low  temperature  end,  this  work  being  about 
2%  lower,  and  excellent  agreement  was  obtained  at  the  highest  tempera- 
tures where  no  difference  was  found.   However  Martin's  results  were 
over  5%  higher  at  the  midpoint,  17°K.  The  differences  along  the  scale 
were  not  random  but  showed  a  clear  divergence  up  to  a  maximum,  and  then 
converged  to  zero  at  the  highest  temperatures.   This  type  of  deviation 
can  most  easily  be  explained  as  some  sort  of  systematic  error  in  the 
determination  of  either  temperature  scale,  or  possibly  a  drift  of  the 
resistance  values  of  the  thermometers.   Martin  discounts  this  latter 
possibility  in  his  work  by  repeated  checks  against  both  vapor  bulb  and 
gas  bulb  thermometery.   In  the  present  work  nearly  every  time  the  sys- 
tem was  at  equilibrium  with  the  helium  bath,  the  thermometers  were 
checked  against  the  vapor  pressure  temperature  of  the  bath.   On  the 
basis  of  this,  it  is  felt  that  there  was  no  resistance  drift  of  the 
thermometers  used  here  at  least  at  the  low  temperatures. 

The  size  of  the  maximum  difference  in  the  two  sets  of  results 
is  quite  large,  but  it  should  be  pointed  out  that  at  17°K  an  error 
in  the  temperature  determination  of  only  1.3%  (.2°)  would  explain 
the  difference.   This  value  of  error  is  outside  of  the  estimated 
error  in  the  temperature  scale  used  here,  however,  and  so  it  was  de- 
cided that  a  recalibration  of  the  thermometers  would  be  performed  at 
the  conclusion  of  the  specific  heat  measurements  on  polyethylene 
(the  copper  data  were  taken  about  half  way  through  the  polyethylene 
work).   It  was  felt  that  with  the  low  temperature  techniques  acquired 
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by  this  experimenter  that  a  more  precise  calibration  might  be  achieved. 
The  recalibration  is  discussed  in  Appendix  A. 

After  the  recalibration  was  performed,  the  heat  capacity  data 
for  copper  were  recalculated,  and  again  compared  with  the  results  of 
Martin.   The  same  relative  deviation  persisted,  however  it  was  approx- 
imately halved  at  the  17°K  point  (2.6%  low).   At  the  highest  temper- 
atures investigated,  30°K,  the  results  were  again  the  same  as  Martin's, 
but  the  same  2%  difference  was  still  present  at  5°K. 

The  new  calibration  was  decided  to  be  the  more  accurate  and  was 
used  based  on  the  following  considerations.   First  the  agreement  with 
Martin's  copper  data  was  better.   Next  the  recalibration  was  per- 
formed after  this  experimenter  had  considerably  more  experience  with 
the  low  temperature  techniques,  and  was  able  to  look  for  small  sys- 
tematic errors  which  may  have  been  overlooked  the  first  time.   Lastly 
this  calibration  was  continuous  over  the  range  2.4  to  30°K  which  was 
highly  desirable  (see  Appendix  A). 

At  the  conclusion  of  this  work,  Alhers  published  specific  heat 
data  on  copper  over  essentially  the  same  temperature  range  as  that 
under  consideration  [ij.   His  work  compares  very  well  with  that  of 
Martin,  but  has  the  same  relatives."  deviation  from  this  work.   For 
comparison,  the  present  copper  specific  heat  data  are  presented  be- 
low along  with  that  of  both  Martin  and  Alhers.   Additionally  the 
parameters  for  the  fitted  curve  equation  (3.2)  are  given,  along  with 
its  o/o  RMS  deviation  from  the  actual  data  points,  (.69%). 
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Parameters   for  the  specific  heat-temperature  fitted  relation 

7  a- I 

C  =  ^  ^a  ~TA  (C    in  mj/°K-gm).      For  T   between   2.4  and    30°K. 


Al 

-.12463X10U 

A2 

.96209X10"1 

A3 

-.23948X10"1 

A4 

.41835X10-2 

A5 

-.2 6 648X10" 3 

A6" 

•  10  385X10-4 

A7 

-.14265X10-6 

%  RMS 

.69 

Deviation 

Molecular  weight  of  copper   63.57  gms/mole 
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COMPARISON  OF  SPECIFIC  HEAT  OF  COPPER 


SPECIFIC  HEATS  IN  MJ/QK/GM 


TEMPERATURE 


C(TUCKEH) 


C(MARTIN) 


C(ALHERS) 


5. 

6. 

1^ 


8. 

9. 

11. 

12. 
JL3^ 
14. 
15. 
16. 
17. 
18. 

19_! 

2o. 
21. 

23. 

24. 
25. 

26. 
27. 
28. 


2V. 

30. 


.144 
.222 

L328 

.465 
.640 

,857 

1.122 

1.443 

1.826 

2.282 

2.820 

3.452 

4.189 

5.043 

_6.028 

7.154 

8.432 

_9_,8  7o 

11.473 

13.243 
15.177 
17.264 
19.489 
21.826 

24.238 
26.679 


.1475 
.2265 
.3334 

.4729 
.6505 

,8718 

1.1444 
1.4737 
1.8708 
2.3409 
2.896o 
3.5440 

4.3oo6 
5.1619 
6.1542 
7.2829 
8.5537 
9.9695 
11.5367 

13.2685 
15.1584 

17.1997 
19.3727 
21.6576 

24.0413 
26.5238 


.1484 
.2275 
.3339 
.4725 
.6488 
.8690 
.1399 
.4690 
.8647 
.3356 
.8910 
.54oo 
.2917 

.1548 
.1388 
.2549 

,5193 


*  INDICATES  NO  DATA 


PI 


APPENDIX   C    -    Raw    Data 


MARLEX     1 


SPECIFIC    HEATS    IN    MJ/OK/GM 
TEMPER  A  TORT  ~XT 

2.597  .18796 


-Z7b~0^  .  1  9 1  54 

2.587  '  .19166 

2.741  .23133 

-2TT6~8 .23560 

2.754  .23643 

2.896  .27448 

2.913 .28187 

2.934  .28523 

3.DU8  .32351 

-3Y0T4 .  33065 

3.098  .33663 

3.200  .37506 


-3723  4 

3.258  .39376 

3.379  .44463 

3.433 .45855 

3.415  .46099 

3.594  .52877 

"3T6~1~8  .52958 

3.617  .53548 

3.808  .61506 

-3Y8tT4 .63022 

3.819  .63708 

3.998  .71218 

TT7003 .72960 

4.046  .75817 

4.196  .82776 

"4T2TO .85138 

4.271  .88733 

4.326  .92379 

~4T340 .93608 

4.366  .96433 

4.419  .98215 

-47  490 1.04819 

4.517  1.05959 

4.539  1.07225 

-4T6  37 T.  1  3612 

4.635  1.14261 

4.673  1.16034 

4T6T9  6 TV  18922 

4.774  1.21758 

4.824  1.29738 

^T&4^ r.  30108 

4.863  1.32905 


•JB- 


MARLEX  1 


SPECIFIC  HEATS  IN  MJ/OK/GM 
TEMPERATURE C 


4.970  1.39622 

-tr^<n : 1 V  4  4  9  5  8 


5.011  "  1.45996 

5. IOC  1.54250 

5.135 r- 1.6 1394 

5.210  1.64436 

5.229  '    1.66338 

5.304 1  .744  34 

5.355  1.79458 

5.418  1.84677 

5.4  68- ! 1  .91434 

5.520  1.96849 

5.604  ,   2.06463 

-SrfrST 2v0  86  88 


5.686  2.16203 

5.746  2.21610 

5.771 2.27986 

5.911  2.42215 

5.930  2.42923 

-6-^)8-6 2.62  558 

5.992  2.62725 

6.181  2.75895 

-6.2  33 2.81262 

6.236  2.84967 

6.513  3.19788 


-6t532 3 .24379 

6.549  3.30149 

6.572  3.31813 

-fr.-Q  6-8 3 .72186 

6.905  3.82183 

6.953  3.85513 

6.892 3.89666 

7.178  4.22355 

7.204  4.23521 

-7 .  22  9 4 .32056 

7.344  4.50184 

7.456  4.63985 

7.535 ■■ 4.83053 


7.541  4.85262 

7.711  4.97886 

-7-.  8  47 5.41  547- 

7.909  5.51952 

7.944  5.53639 


-8T-V57 5^94872 

8.208  6.12060 


93 


MARLEX  1 


SPECIFIC  HEATS  IN  MJ/OK/GM 


"mrMFERTATUSnE" 


8.277   6.29253 

"87437  6.53970 


8.671  *  7.00062 

8.623  7.19002 

8.773 7.29157 

8.985  7.70739 

9.075  7.98402 

-97077 8.03411 

9.353  6.54090 

9.415  8.88594 

"97424 "8793725 

9.610  9.34961 

9.820  9.74374 


~ 9T81+7 9.97924 

,    9.990  10.618U3 

10.265  11.22138 

10.331 11.48649 

10.604  12.30049 

10.834  12.99738 

T0T892 15.1  7960 

11.148  14.03310 

11.339  14.67437 


7774 CTO 14.93354 

11.644  15.82448 

11.786  16.29977 

11,949 —  16.91263 

12.162  17.62685 

12.198  17.79657 

727832 2  0728792 

12.859  20.48844 

13.239  22.24001 

1  376  42 2  3.72455 

13.635  23.74350 

13.856  24.79318 


7  47326 2  6.70627 

14.343  26.93295 

14.563  27.92274 

~147T988 29.89960 

14.973  31.75482 

15.371  31.86838 

T57634 32.10680 

16.320  37.06861 

17.374  42.64875 

7774~83 43.70005 

17.819  45.38054 
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MARLEX     1 


SPECIFIC    HEATS    IN    MJ/OK/GM 
-TEMPER7UUTE ~X~ 

18.430  48.82613 

r8v4  6  6 49.36476 

18.875  '  51.37767 

19.274  54.50185 

19.555 "56.01417 

19.799  57.71516 

20.105  59.66975 

20.721* 63.97106 

20.884  65.13815 

21.133  68.59212 

21.667 71.31811- 

21.910  71.80872 

22.468  77.27470 

23.028 8n.64122 


22.828  81.65808 

23.U99  85.04005 

23.845 88.50399 

24.186  91.55677 

24.980  95.85445 

-25.00  1 98.  1  4700 

25.292  100.16850 

26.191  106.23490 

26.578 115.46438 

27.159  118.91252 

27.572  121.90134 

^8r8^5 132.66855 
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MARLEX  2 


SPECIFIC  HEATS  IN  MJ/OK/GM 


TEMPER  ATO^T 


4.349 

"4.375 
4.373 
4.538 
4.540 
4.568 
4.703 

"4771*0" 
4.762 
4.881 

TJT94"H" 
4.967 
5.092 

"57170" 
5.187 
5.301 

"57390- 
5.484 
5.412 

-5T6^T6" 
5.670 
5.645 

"57850" 


5.859 
5.897 

"67078" 
6.197 
6.165 

-67378" 
6.454 
6.595 

-67TT8- 
6.770 
6.952 

77030" 
7.118 
7.315 

"773  30 
7.489 
7.642 

T76TO" 
7.875 
7.990 

-8r098~ 
8.338 


.10836 
.15563 
.15943 
.26999 
.28241 
.31700 
.41884 
.46093 
.50724 
.60343 
.66051 
.70458 
.82800 
.91897 
.96224 
.05350 
.17302 
.22102 
.25201 
.44090 
.51030 
.53578 
.77896 
.78156 
.89658 
.09759 
.28717 
.31456 
.60245 
.75170 
.93908 
.14388 
.29967 
.57285 
.73857 
.93654 
.26753 
.29696 
.69568 
.90489 
.04339 
.53743 
.69330 
.94474 
.47455 
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MARLEX  2 


SPECIFIC  HEATS  IN  MJ/OK/GM 
TEMPER  ATU3TE C 


8.612  8.36381 

— 8.6  73 8.37944 

8.720  '  8.42132 

9.164  9.63058 

— 9.337 -  10.13203 

9.447  10.62664 

9.711  11.61380 

—9.907 1  1 .85851 

10.000  12.28690 

10.314  13.08954 
-10.4  56 13.67324 

10.534  14.10426 

10.860  15.03178 

-11.017 15.61562 

11.098  16.19879 

11.391  16.97032 

-1  1 .600 1  7.  80 1  59 

11.694  18.33617 

11.938  19.13534 

-12.195 20.23688 

12.315  20.92201 
12.493  21.41557 

—12.797 22.72404 

12.961  23.71420 

13.066  23.89595 

-13.4  1  3 25.4327  1 

13.701  26.84314 

13.660  26.95160 

-VU--.  144 29.02634 

14.353  30.01182 

14.M22  30.75287 

-14.929 32.89604 

15.030  33.22668 

15.236  34.73176 


-1-5-.-T3-8" 3  7.14855 

15.784  37.24379 

16.060  39.29161 

1  fr.  5  73 41.48706 

16.596  42.02919 

16.911  44.06952 

17.397 4  6  .  2492  8 

17.448  46.91189 

17.827  49.84552 

1 8  .  1  45 5  0.75218 

18.345  52.51535 
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MARLEX  2 


SPECIFIC  HEATS  IN  MJ/OK/GM 

TEMPERATURE  "TT" 

18.832  55.59494 

" 18. /89  55.66262 

19.311  '  58.92028 

19.688  62.097U1 

19.T9T) 62.71279 

20.32U  66.22665 

20.554  68.59343 

20 .833 71.12623 

21.507  73.83693 

21.511  75.11742 

21.938 78.73929 


22.729  82.89937 

22.647  84.30054 


-2370875 88 .  27922 

23.871  92.26345 

23.844  96.20608 

"2  472"56 96.79725 

24.965  102.41275 

25.217  104.56977 


-257*6-9 110.32250 

26.118  112.90491 

26.716  117.49972 

-2-67707 1 1  9.71567 

27.999  131.33693 

28.399  131.63365 


"257-6171 1  32.40289 

29.244  137.68334 
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LOW  DENSITY 


SPECIFIC  HEATS  IN  MJ/OK/GM 


TEMPERATURE 


2.521 
"2.533 

2.579 
2.621 

"2.655- 
2.695 
2.731 

-2t7^0- 
2.816 
2.851 

-2v9tT8- 
2.9U1 
2.978 
3.063 
3.097 
3.088 
3.241 
3.251 
3.266 
3.454 
3.U49 
3.440 

-3.645 
3.6U8 
3.657 
3.843 
3.851 
3.863 

-*nrt>^2- 
4.061 
4.064 
4.278 
4.328 
4.318 
4.335 
4.U74 
4.498 
4.497 
4.508 
4.642 

HT.-66-2- 


.28781 
V29686 
.31867 
.33245 
.35445 
.37404 
.38825 
.41757- 
.42907 
.44835 
.47711 
.50078 
.51675 
.57180 
.57700 
.59378 
.65976 
.67721 
.69384 
.81479 
.81595 
.82114 
.96564 
.96623 
.97461 
13362 
14367 
14469 
33429 
33820 


4.706 

4.831 

-4^888" 

4.932 


1 .34516 
1.57952 
1.61775 
1.62234 
1.63967 
1.81021 
1.81792 

■1.84558- 
1.86155 
2.03621 

-2vO6809^ 
2.13112 
2.30063 

-2x38673 
2.46734 
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LOW    DENSITY 


SPECIFIC    HEATS    IN    MJ/OK/GM 


T^PEPTATTOE  ~X 

5.034  2.63293 


37TT6  2T  7  4  9  0  4 

5.140  '  2.79630 

5.189  2.88276 

5.245 2.9665  9 

5.360  3.19051 

5.382  3.21864 

T74~5o 3.36224 


5.470 
5.615 

3.36955 
3.64384  • 

5.634 
5.712 
5.713 

3.68241 
3.83U01 
3.85035 

5.895 
;                                              5.916 

5.980 

4.22231 
4.27282 
4.32430 

5.99Z 
6.209 
6.235 

4.45900 
4.91476 
4.95730 

6.303 
6.550 

5.03374 
5.10324 
5.69680 

6.5/5 
6.592 
6.644 

5.76155 
5.78959 
6.00095 

6.904 
6.928 
6.922 

6*59770 

6.63617 
6.65371 

7.268 

"7.52172 

7.263  7.56141 

7.287  7.61844 

7.625 8  .  57451 

7.669  8.71615 

8.024  9.80153 

~8 .  08"2 9.9  8  70  T 

8.405  11.06141 

8.460  11.21426 


-8T5D5 1 1  .  36980^ 

8.815  12.46119 

8.899  12.75976 

"8T933 1 2.86  68  7~ 

9.231  13.96903 

9.361  14.44571 

"9T3B7 T4.58159 

9.661  15.61894 
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LOW    DEMSITY 


SPECIFIC    HEATS    IN    MJ/OK/GM 

TEMPER  ATWE ~C~ 

9,884  16.52142 

9Y894 16.  57931 

10.127  17.59983 

10.464  18.916U2 

10.452 18.95116 

10.639  19.73098 

11.038  21.57511 

11.064 21  .69968 

11.194  22.29614 

11.631  24. 39704 


-1-1-V6  69 2  4.42918 

11.788  25.18032 

12.229  27.30692 

T2r27ib 27.65680 

12.410  28.33392 

12.845  30.68352 

-1  2 .90  1 30.  86788 

13.063  51.55914 

13.515  34.28889 

-1-3T5T6 34.61570 

13.767  35.76025 

14.245  38.38456 

14.31  1 3  8.66656 

14.528  39.96585 

15.022  42.96486 

T^5v094 43.32534 

15.340  44.96633 

15.862  48.09940 

15.940 48.4751  1 

16.220  50.30252 

16.758  53.82994 

16.845 54.36106 

17.137  56.28382 

17.700  59.96677 

17.785 60.70527 

18.713  67.33329 

18.777  67.69524 

19.002 69.52506 

19.786  75.32322 

19.360  76.15117 

19.986 76.79421- 

20.904  84.22005 

20.958  84.50975 

-2-K-fr3* 85-.-0S1 82 

22.047  92.92935 
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LOW    DENSITY 


SPECIFIC    HEATS    IN    MJ/OK/GM 
TEWETTATWE  ~C~~ 

22.118        94.18537 

*:.     22.118  94.42999 

23.228  103.33531 

23.336  103.U3852 

23.211 "103.96273 

24.385  113.53377 

24.352  115.23389 

■ — -2UVS6T) 1 16.08967 

25.490  124.87124 

25.557  125.00688 

25".7T6 125.96857 

26.735  133.83623 

26.65U  134.54346 

T6V992 135.14483 


27.869  140.81310 

27.939  143.93052 

"Z872T5 1 4  5.08745 

29.122  154.61557 

29.283  156.00964 


29.038 161.87336 
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THE     100    0/0    CRYSTALINE    SAMPLE 


HEAT    CAPACITIES    IN    MJ/OK/GM 


TEMPERATURE 


2, 

,5 

3, 

,  0 

3, 

,5 

4, 

,0 

4. 

.  5 

5. 

,  0 

5, 

,  5 

6, 

,0 

6. 

,5 

7, 

,0 

7, 

,5 

8. 

,0 

8, 

,5 

9, 

,0 

9, 

,  b 

10. 

,0 

10. 

,  5 

11. 

,0 

11. 

,  5 

12, 

,0 

12. 

5 

13. 

0 

13. 

5 

14, 

,0 

14. 

5 

15. 

0 

15. 

,5 

16, 

,0 

16, 

,5 

17. 

0 

17. 

5 

18, 

,0 

18. 

,5 

19. 

,0 

19. 

5 

20. 

,G 

20. 

5 

21. 

0 

21. 

5 

22. 

0 

22. 

5 

23. 

0 

23. 

c 

24. 

0 

24. 

5 

25. 

0 

25. 

5 

26. 

0 

26. 

5 

27. 

0 

50. 

,0 

100. 

0 

150. 

0 

200. 

0 

C  ACTUAL 

TARASOV  C 

PER  CENT 
DIFFERENCE 

.124 

.123 

.9 

.211 

.213 

-.7 

.335 

.337 

-.9 

.501 

.504 

-.6 

.715 

.717 

-.3 

.984 

.984 

.0 

1.31? 

1.310 

.3 

1  .708 

1.700 

.5 

2.174 

2.162 

.5 

2.715 

2.700 

.6 

3.337 

3.321 

.5 

4.045 

4.0  50 

.4 

4.842 

4.832 

.2 

5.732 

5.735 

-.1 

6.719 

6.742 

-.3 

7.805 

7.858 

-.7 

8.994 

9.i)87 

-1.0 

10.288 

1  ?.433 

-1  .4 

1  1.689 

1 1.898 

-1.8 

13.197  , 

13.485 

-2.2 

14.815 

15.195 

-2.6 

16.543 

17.029 

-2.9 

18.381 

18.985 

~Z,3 

20.330 

21 .064 

-3.6 

22.389 

2  3.26  3 

-3.9 

•24.558 

25.581 

-4.2 

26.836 

28.013 

-4.4 

29.221 

30.557 

-4.6 

31.71?   % 

33.210 

-4.7 

34.310 

35.966 

-4.8 

37.010 

38.822 

-4.9 

39.812 

41.773 

-4.9 

42.714 

44.815 

-4.9 

45.715 

47.942 

-4.9 

43.813 

51.151 

-4.8 

52.007 

54.435 

-4.7 

55.295 

57.792 

-4.5 

58.677 

61.215 

-4.3 

62.152 

64.701 

-4.1 

65.720 

68.245 

-3.8 

69.381 

71 .844 

-3.5 

73.136 

75.492 

-3.2 

76.990 

79.187 

-2.9 

■80.942 

82.924 

-2.4 

84.996 

86.701 

-2.0 

89.156 

90.513 

-1.5 

9^.429 

94.358 

-1.0 

97.821 

98.233 

-.4 

102.339 

102.135 

.2 

106.994 

106.061 

.9 

332.800 

290.367 

12.8 

682.400 

644.550 

5.5 

942.900 

920.149 

2.4 

1197.000 

1117.313 

6.7 
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THE    MARLEX    ONE    SAMPLE 


HEAT    CAPACITIES    IN    MJ/OK/GM 


TEMPERATURE 


5 
0 
5 
0 

5 
0 
5 
0 
5 


2, 

3. 

3. 

4. 

4, 

5, 

5. 

6. 

6. 

7.0 

7.5 

8.0 

8.5 

9.0 

9.5 

10.0 

10.5 

11.0 

11.5 

12.0 

12.5 

13.0 

13.5 

14.0 

14.5 

15.0 

15.5 

16.0 

16.5 

17.0 

17.5 

18.0 

18.5 

19.0 

19.5 

20 

20. 

21 

21 

22 

22. 

23.0 

23.5 

24.0 

24.5 

25.0 

25.5 

26.0 

26.5 

27.0 

50.0 

100.0 

150.0 

200.0 


0 
5 
0 
5 
0 
5 


C  ACTUAL 

TARASOV  C 

PER  CENT 
DIFFERENCE 

.172 

.171 

.5 

.303 

.303 

-•9 

.487   , 

.490 

-  .5 

.735 

.759 

-.5 

1  .05* 

1  .057 

-.3 

1  .449 

1.448 

.1 

1.928 

1.920 

.4 

2.495 

2.478 

.7 

•  3.154 

3.127 

.8 

3.909 

3.875 

.9 

4.761 

4.725 

.8 

5.714 

5.682 

.6 

6.76  3 

6.749 

.5 

7.926 

7.932 

-.1 

9.  187 

9.231 

-.5 

10.553 

10.64  7 

-.9 

12.u2  3 

12.184 

-1.3 

13.597 

13.84  1 

-1.8 

15.275 

15.618 

-2.2 

17.057 

17.515 

-2.7 

18.94  3 

19.530 

-3.] 

20.931 

21  .663 

-3.5 

2  3.021 

23.911 

-5.9 

25.213 

26.272 

-4.2 

27.506 

28.742 

-4.5 

29.899 

31.319 

-4.7 

32.392 

35.998 

-5.0 

34.985 

36.776 

-5.1 

37.676 

39.648 

-5.2 

40.465 

42.611 

-5.3 

45.353 

45.659 

-5.3 

46.338 

48.789 

-5.3 

4  9 .  4  2  0 

51.995 

-5.2 

52.600 

55.274 

-5.  i 

55.876 

55.621 

-4.9 

59.248 

62.031 

-4.7 

62.717 

65.501 

-4.4 

66.281 

69.026 

-4.] 

69.940 

72.602 

-3.8 

73.693 

76.226 

-3.4 

77.539 

79.893 

-3.0 

81.476 

83.601 

-2.6 

85.505 

87.346 

-2.2 

89.622 

91 .125 

-1.7 

93.825 

94.934 

-1.2 

98.  112 

98.772 

-.7 

102.479 

102.634 

-.2 

106.92  3 

106.520 

.4 

11 1.439 

110.426 

.9 

1 16.021 

114.350 

1.4 

332.800 

296.206 

11.0 

682.400 

646.873 

5.2 

942.900 

921  .362 

2.3 

1 197.000 

1118.064 

6.6 
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THE    MARLEX    TWO    SAMPLE 


HEAT 

CAPACITIES 

IN  MJ/OK/GM 

EMPERATURE 

C  ACTUAL 

TARASOV  C 

PER  CENT 
DIFFERENCE 

2.5 

» 

.201 

• 

3.0 

« 

.357 

• 

3.5 

* 

.582 

« 

4.0 

* 

.831 

* 

4.5 

1  .249 

1.261 

-1.0 

5.0 

1.735 

1  .729 

.4 

5.5 

2.315 

2.288 

1.2 

6.0 

2.995 

2.947 

1.6 

6.5 

3.778 

3.711 

1.3 

7.0 

4.666 

4.584 

1.8 

7.5 

5.66? 

5.573 

1.6 

8.0 

6.765 

6.679 

1.3 

8.5 

7.977 

7.907 

.9 

9.0 

9.297 

9.258 

.4 

9.5 

10.724 

10.733 

-.1 

10.0 

12.258 

12.332 

-.6 

10.5 

13.898 

14.054 

-1.1 

11.0 

15.643 

15.398 

-1  .6 

11.5 

17.490 

17.863 

-2.1 

12.0 

19.440 

19.947 

-2.6 

12.5 

21.491 

22.148 

-3.1 

13.0 

23.642 

24.462 

-5.5 

13.5 

25.891 

26.886 

-3.8 

14.0 

2  8.24  0 

29.417 

-4.2 

14.5 

3  J. 685 

32.051 

-4.4 

15.0 

33.229 

34.784 

-4.7 

15.5 

35.869 

37.612 

-4.9 

16.0 

38.606 

40.531 

-5.0 

16.5 

41.439 

43.536 

-5.1 

17.0 

44.370 

46.623 

-5.1 

17.5 

47.598 

49.787 

-5.0 

18.0 

50.522 

53.024 

-5.0 

18.5 

53.744 

56.330 

-4.0 

19.0 

57.062 

59.701 

-4.6 

19.5 

60.476 

63.131 

-4.4 

20.0 

63.985 

66.617 

-4.1 

20.5 

67.588 

70.155 

-3.8 

21.0 

71.283 

73.742 

-3.4 

21.5 

7  5.0  67 

77.372 

-3.1 

22.0 

78.956 

81 .044 

-2.7 

22.5 

82.887 

84.753 

-2.5 

23.0 

86.913 

88.497 

-1.8 

23.5 

91 .007 

92.272 

-1.4 

24.0 

95.161 

96.076 

-1.0 

24.5 

99.365 

99.905 

-.5 

25.0 

105.607 

103.758 

-.1 

25.5 

107.872 

107.631 

.2 

26.0 

112.144 

11 1 .523 

.6 

26.5 

116.405 

115.432 

.8 

27.0 

120.633 

119.355 

1.1 

50.0 

332.800 

299.732 

9.9 

100.0 

682.400 

648.275 

5.0 

150.0 

942.900 

922.095 

2.2 

200.0 

1197.000 

1118.518 

6.6 
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THE    LOW    DENSITY    SAMPLE 


HEAT    CAPACITIES    IN    MJ/OK/GM 


TEMPERATURE 


2. 

,5 

3. 

0 

3. 

5 

4. 

0 

!»« 

,5 

5. 

0 

5. 

b 

6. 

0 

6. 

,5 

7. 

0 

7, 

5 

8. 

0 

8, 

,5 

9, 

0 

9. 

,5 

10. 

,0 

10, 

,5 

11, 

,0 

11, 

5 

12. 

,'J 

12. 

5 

13, 

0 

13, 

,5 

14, 

,0 

14. 

,  5 

15, 

0 

15. 

,  5 

16. 

,  0 

16. 

5 

17. 

,0 

17. 

,5 

18. 

,0 

18, 

,5 

19, 

,0 

19, 

.5 

20. 

,0 

20. 

,5 

21, 

,  0 

21, 

,  5 

22, 

.0 

22. 

,  5 

23. 

,0 

23. 

,5 

24, 

.0 

24, 

,5 

25, 

,0 

25, 

.5 

26, 

,0 

26, 

,5 

27, 

,  0 

50, 

,0 

100, 

,  0 

150, 

,0 

200, 

,0 

C  ACTUAL 

TARASOV  C 

PER  CENT 
DIFFERENCE 

.288 

.288 

.0 

.523 

.520 

.6 

.855 

.855 

-.1 

1.298 

1.304 

-.5 

1  .865 

1  .871 

-.3 

2.566 

2.564 

•  i 

3.405 

3.386 

.6 

4.386 

4.345 

.9 

5.510 

5.448 

1.1 

6.776 

6.697 

1.2 

8.  182 

8.098 

1.0 

9.724 

9.651 

.7 

11 .397 

1 1 .356 

.4 

13.197 

13,210 

-.1 

13.118 

15.209 

-.6 

17.15? 

17.349 

-1.1 

19.293 

19.624 

-1.7 

21.546 

22.028 

-2.2 

2  3.892 

24.554 

-2.8 

26.331 

27.195 

-3.3 

28.86^ 

29.945 

-3.0 

31.473 

32.798 

-4.2 

34.169 

35.746 

-4.6 

36.9U5 

38.784 

-5.0 

39.799 

41 .906 

-5.3 

42.732 

45. 106 

-5.6 

45.742 

48.578 

-5.8 

U8.832 

51 .716 

-5.9 

5  2.002 

55.  117 

-6.0 

55.255 

58.575 

-6.0 

58.592 

62.085 

-6.0 

62.017 

65.643 

-5.8 

65.532 

69.246 

-5.7 

69. 141 

72.888 

-5.4 

72.845 

76.568 

-5.1 

76.647 

80.280 

-4.7 

80.549 

34.022 

-4.3 

84.551 

87.791 

-3.8 

88.651 

91.584 

-3.3 

92.848 

95.399 

-2.7 

97. 137 

99.233 

-2.2 

101.511 

105.083 

-1.5 

105.961 

106.948 

-.9 

1 10.476 

1U-. 826 

-.3 

115.038 

1  14.714 

.3 

119.628 

118.612 

.8 

124.222 

122.517 

1.4 

128.791 

126.429 

1.8 

133.300 

150.345 

2.2 

137.710 

134.265 

2.5 

332.800 

310.235 

6.3 

682.400 

652.453 

4.4 

942.900 

924.277 

2.0 

1197.000 

1119.869 

6.4 
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THE  0  0/0  CRYSTALINE  (AMORPHOUS)  SAMPLE 


HEAT  CAPACITIES  IN  MJ/0K/3M 


TEMPERATURE 


2, 

.  5 

3, 

,0 

3. 

,5 

U. 

,0 

4, 

.  5 

5, 

,0 

5, 

■  5 

6. 

,0 

6. 

,5 

7. 

,0 

7. 

,  5 

6. 

,0 

8. 

,5 

9. 

,  0 

9, 

5 

10, 

,0 

10. 

,  5 

11. 

0 

11. 

5 

12. 

,0 

12. 

5 

13. 

0 

13. 

,5 

14. 

,0 

14. 

5 

15. 

0 

15. 

5 

16. 

.0 

16. 

,5 

17. 

0 

17. 

5 

18. 

,0 

18. 

5 

19. 

0 

19. 

5 

20. 

0 

20. 

5 

21. 

0 

21. 

,5 

22. 

0 

22. 

5 

23. 

0 

23. 

5 

24. 

0 

24. 

5 

25. 

0 

25, 

5 

26. 

0 

26. 

5 

27. 

0 

50. 

0 

100. 

0 

150. 

0 

200. 

0 

C  ACTUAL 

TARASOV  C 

PER  CENT 
DIFFERENCE 

.427 

.428 

_  n 

.78P 

.781 

1.0 

1.296 

1.295 

.1 

1.97U 

1.983 

-.4 

2.841 

2.851 

-.3 

3.90P 

3.904 

.1 

5.179 

5.148 

.6 

6.658 

6.590 

1.0 

8.341 

8.235 

1.3 

10.222 

10.039 

1.3 

12.292 

12.152 

1.1 

14*542 

14.421 

.8 

16.959 

16.891 

.4 

19.531 

19.552 

-.1 

22.244 

22.394 

-.7 

25.084 

25.403 

-1.3 

28.040 

28.565 

-1.9 

31.097 

31.866 

-2.5 

34.246 

35.291 

-3.1 

37.475 

58.827 

-3.6 

40.776 

4  2.460 

-4.1 

44.140 

46.177 

-4.6 

47.563 

49.967 

-5.1 

•  51 .041 

53.819 

-5.4 

54.570 

57.723 

-5.8 

58.151 

61.672 

-6.1 

61.783 

65.655 

-6.3 

65.471 

69.668 

-6.4 

69.216 

73.704 

-6.5 

73.025 

77.756 

-6.t> 

76.903 

81.822 

-6.4 

80.856 

85.895 

-6.2 

84.891 

89.974 

-6.0 

89.015 

94.053 

-5.7 

95.234 

98.132 

-5.3 

97.553 

102.207 

-4.8 

101.975 

106.277 

-4.2 

106.502 

110.339 

-3.6 

111.13? 

114.593 

-2.9 

115.864 

118.437 

-2.2 

120.685 

122.470 

-1.5 

125.584 

126.492 

-.7 

13J.541 

130.501 

.0 

135.533 

134.498 

.8 

140.527 

138.481 

1  .5 

145.481 

142.451 

2.1 

150.34  8 

146.408 

2.6 

155.068 

150.351 

3.0 

159.569 

154.280 

3.3 

163.770 

158.195 

3.4 

332.800 

327.092 

1.7 

682.400 

659.158 

3.4 

942.900 

927.779 

1.6 

1 197.000 

1122.037 

6.3 
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